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Abstract 

The SU(5) unification-scale relation hj, = h T between the 6-quark and r- 
lepton Yukawa couplings severely constrains tan (5 and the f-quark mass (even 
more so if ht = hb = h T holds) in supersymmetric models. We examine the 
implications of these constraints for the Higgs sector assuming universal soft 
breaking terms, and emphasize that both of these relations impose unique 
characteristics in terms of symmetries and of the spectrum. We further study 
the tan/3 « 1 scenario, which is suggested by hb = h T , and, in particular, the 
loop-induced mass of the light Higgs boson. We compare the effective poten- 
tial and renormalization group methods and stress the two-loop ambiguities 
in the calculation of the mass. These and a large enhancement to the loop cor- 
rection due to t-scalar left-right mixing considerably weaken the upper bound 
on the mass of the light Higgs boson that has been reported. Nevertheless, we 
find that for this scenario the Higgs boson is probably lighter than 110 GeV, 
and typically lighter than 100 GeV. Thus, it is in the mass range that may 
be relevant for LEPII. Our numerical results are presented in a self-contained 
manner in section In separate appendices we discuss the global symmetries 
of the Higgs potential, the issue of false (color-breaking) vacua, which may be 
important for tan/3 ~ 1, two-loop calculations, and the effect of an additional 
Higgs singlet. We show that the approximate constraints that are often used 
to eliminate color-breaking vacua are not always relevant. 

PACS numbers: 12.10.Dm, 11.30.Pb, 14.80.Gt, 14.80.Ly 
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I. INTRODUCTION 



Supersymmetric grand- unified theories JTJJ^] (SGUT), though far from being proved, are 
an attractive framework on both observational and theoretical grounds. (Their relation 
to string theory remains at present poorly understood.) In the minimal models, which 
we assume hereafter, there is a grand desert between the weak and the unification scales 
(M z < Q < TeV and M G « 10 16 - 10 17 GeV, respectively). At the unification point M G 

oci(Mq) = a 2 (M G ) = a 3 (M G ) = a G (1) 

(where a\ = |ay and a G ~ 0.04), and, depending on the grand unifying group representa- 
tions of the Higgs superfields that couple to matter, certain relations between the standard 
model Yukawa couplings hold. In the minimal models based on SU(5), SO(10), and E 6 
these representations are the fundamental ones, and one obtains || 

h T (M G ) = h b (M G ). (2) 

For some models based on SO (10) (or larger groups) one has the stronger prediction 

h T (M G ) = h b (M G ) = ht(M G ), (3) 

but this only holds if one makes the additional assumption that the masses are generated by a 
single complex Higgs 10-plet. (h T b t are the SM Yukawa couplings of the r, b, t-superfields, 
respectively.) One usually assumes that some perturbation modifies the coupling or the 
masses of the two light families where, in principle, similar relations should, but do not, 
holdQ. (See also Appendix |D|.) Also, such relations often do not hold in superstring unified 
models. 

At the weak-scale the particle content is that of the minimal supersymmetric standard 
model (MSSM), i.e., that of a two Higgs doublet model but with each bosonic (fermionic) 
degree of freedom complemented by a fermionic (bosonic) one. In particular, the Higgs 
sector contains three Goldstone bosons and five physical degrees of freedom - two CP even 
(h° and H°) and one CP odd (A ) neutral, and one complex charged (H + ) Higgs bosons. 
The symmetries and spectrum of the Higgs scalars are strongly affected by the above GUT 
assumptions, and are the subject of our present study. However, a review of previous results 
in SGUT's, which establish the base for our present work, is appropriate. 

The consistency of ([!]) with the data was pointed out recently in Ref. [Q and that of 
(0) and @ in Ref. ||. Relation ([]]), however, holds only if there is exactly one pair of 
Higgs doublets (but any number of SM singlets is allowed). We previously investigated the 
status of coupling constant unification in great detail || where we considered non-trivial 
matching conditions^ that perturb relation (|l|) in a model dependent way. We concluded 



1 Since the Yukawa couplings corresponding to the third family are much larger than those of the 
two light families, relation (^) [or (||)] will be only slightly affected. 

2 Heavy degrees of freedom (which are model dependent) are integrated out of the effective theory 
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that no significant constraints on the MSSM parameters can be deduced just by that relation. 
However, we findP] that ([!]) predicts for the strong coupling [using a{Mz)~ x = 127.9 ±0.1 
and s 2 {M z ) = 0.2324 ± 0.0003 (for m p t ole = 143 GeV) for the (MS) fine-structure constant 
and weak angle, respectively] 

a s (M z ) = 0.125 ± 0.001 ± 0.008 + H as + 3.2 x 10- 7 GeV~ 2 [(mf ie ) 2 - (143 GeV) 2 ] , (4) 

where the ±0.008 (±0.001) uncertainty in (|4]) is theoretical (due to the input parameter error 
bars). The former is due to the unknown values of the matching conditions (i.e., threshold 
and nonrenormalizable-operator effects) mentioned above. The function H as is a correction 
from Yukawa interactions. It is negative but negligible, unless some Yukawa couplings are 
> 1, for which we obtain H ag « —(0.001 — 0.003). The quadratic m pole dependence is 
induced by the correlation between the weak angle extracted from the data and the t-quark 
pole mass, m pole . Thus, typically one predicts a s (M z ) > 0.12 with a significant dependence 
on m pole . This observation is of importance when simultaneously considering ([I]) and (0) [or 
(0)] to predict the 6-quark (MS running) mass m h [and/or the t-quark mass, if using (|3|)] in 
terms of the r-lepton mass. We find [using, e.g., Eq. (||) but with a more careful treatment 
of theoretical uncertainties to the 6-quark mass prediction] that only a small region in the 
tan/5 — m pole plane is allowed by the experimental range of m& ||, where tan/5 = fup/^down 
is the ratio of the two Higgs doublet expectation values. This holds even when non-trivial 
matching conditions at both scales are considered. 

It is customary to assume that the MSSM spectrum is given to a good approximation 
near Mq by a small number of universal supersymmetry-breaking soft parameters; i.e., a 
common scalar mass tuq, trilinear and bilinear (dimension-one) couplings Aq and Bq, and a 
common gaugino mass Mi; in addition to a supersymmetric Higgs mass (i.e., the Higgsino 
mass), fiQ [Q]. We will adopt these assumptions below. (For an alternative scenario, see, 
for example, Ref. §J.) However, our conclusions so far are largely (i.e., except for the the 
details of the matching conditions) independent of those assumptions and hold in any SGUT 
scenario with an approximate grand-desert in which the MSSM /^-functions are valid. The 
universality assumption just strengthens the constraints on the tan/? — m pole plane; i.e., 
assuming universal initial conditions and requiring that the weak-scale MSSM (one-loop 
improved) scalar potential is consistent with a broken SU(2) x U(l) symmetry, we have to 
further constrain tan (3 > 1 and ht > h b . We then find two allowed regions, 

1. tan/3 = 1 + Ap, 140 < m p t ole (< 200) GeV, h t « 1 > h b ; 

2. tan/? « 50 ± 10, 170 ± 10 < mf le (< 200) GeV, h t >h b ^ 1. 



at Mq and some modification to the naive boundary conditions in ([![) and (||) [or (^)] is required, in 
general. The appropriate matching functions may also include corrections from nonrenormalizable 
operators since Mq is not too far from the Planck scale. Similarly, integrating out the MSSM new 
particle spectrum will correct the boundary conditions at M Z - 

3 Eq. (|j) updates the results in || for the effect of more recent precision electroweak data. 



3 



The first of these solutions is consistent only with hb = h T , and has A^ < 1 for m p t ole < 190 
GeV. A part of the region which is described by the second solution is consistent with 
the stronger prediction h T = = h t . A third solution with 3 < tan/3 < 40 exists for 
rrlF t ole ~ 215 ±10 GeV, which is the upper bound on m? t ole if the model is to stay perturbative 
up to Mq. However, such large values for m t are inconsistent with precision electroweak 
data. The allowed regions are illustrated in Fig. |l| 

These results are in good agreement with those of Ref. with any small differences 

due to different treatments of matching conditions and of the coupling constant unification 
condition. Arbitrarily relaxing both condition (]!]) and (Q) independently [e.g., a s (Mz) = 
0.11 and hb(Mc) = 0.85h T (Mc)] could invalidate our conclusions. However, the corrections 
to ([]]) and (0) are strongly correlated, and when treated as such the strong constraints 
on the tan/3 — m? t ole plane holdQ. In particular, the m&/m r ratio is relatively insensitive 
to the details of the decoupling of the low-scale spectrum [|J, and our conclusions would 
not be significantly modified if those thresholds were treated differently than in our (semi- 
analytic) approximation (where we used six effective mass parameters). Let us stress that 
any treatment of theoretical uncertainties would be lacking unless the potentially important 
corrections at the high scale are accounted for in one form or another. These issues are 



summarized in Ref. flqjlq| , where the second reference updates some of the results of Ref. 
HH. Here, relevant results are again updated using a preliminary analysis of most recent 
precision data. 

A two-parameter fit [allowing the light (SM-like) Higgs boson mass, m h o, to vary from 
50 — 150 GeV with a central value m h o = Mz] to all Z, W, and neutral-current data yields 

s 2 (M z ) = 0.2324 ± 0.0003 - 0.92 x 10~ 7 GeV- 2 [(mf' e ) 2 - (143 GeV) 2 ], (5) 

m pole = M3 +17 + U 5 ln GeV; (6) 

which in turn implies tan/3 < 1.3 (at one s.d.). Eqs. (|5|) and (ffl) update Ref. [|H]] where 
rnF° le = 134±28 + 12.5 ln GeV was given (which also implied tan/3 < 1.3). 

It was pointed out by Barger et al. [O that if indeed mf ie < 160 GeV (as suggested by 
the range given in Ref. |17j and here) then m^o [in case (1)] is determined by the magnitude 



of the loop corrections, A/jO, and is ^ 85 GeV. If so, the 6-quark mass prediction in minimal 
SU(5) (or similar) SGUT constrains the SM Higgs boson mass more significantly than the 



general MSSM (triviality) upper bound of ~ 130 GeV [|18H20|1. We agree qualitatively with 
this conclusion, but argue that two-loop ambiguities and important left-right t-scalar mixing 
effects, in addition to the higher mf°' 6 range, weaken the bound to m h o ^110 GeV. 

Below, we will establish the basis for and then pursue further numerical studies of the 
prediction for m^o and its upper bound in the tan /3 w 1 scenario. Rather than approximating 
the one-loop correction we will calculate m^o for a given point in the parameter space using 
the supersymmetric spectrum (calculated numerically), and use monte-carlo routines to 
study the upper bound. In section [TT] we briefly review the MSSM Higgs potential and its 



J That is unless one assumes a conspiracy of various corrections that are otherwise independent. 
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minimization conditions. We specialize the discussion to the minimal SGUT scenarios, cases 
(1) and (2), and emphasize that the Higgs potential exhibits approximate global symmetries 
in these cases, on which we elaborate in Appendix [A]. We concentrate thereafter on the 



tan (3 « 1 scenario. Some general features are described in section [TTT[ In particular, the 



large Higgsino mass parameter can generate false vacuum solutions with broken color and 
charge. We discuss that issue in greater detail in Appendix ||. The calculation of the 
light Higgs boson mass in both the effective potential and renormalization group methods 
is discussed in section |IV|, where we use the comparison between the two methods to study 
the two-loop ambiguities in the calculation. The details of the two methods, as well as 
two-loop calculations, are further discussed in Appendix Section [V] is reserved for further 
numerical studies of the Higgs boson mass, and we show that because of two-loop ambiguities 
and ii — tn mixing m^o can easily exceed 85 GeV but, unless m? t ole is much heavier than 
the range suggested by (^|), it remains in the range that may be relevant for LEPII. There 
is even some possibility that it is in the range still relevant for LEPI. We also discuss the 
relation between t-scalar mixing and color breaking. In particular, parameter ranges which 
correspond to large mixing (and large m^o) often have an unacceptable color-breaking global 
minimum, but not always. The discussion in section is to a great extent independent of 
the other sections and can be read on its own. We will summarize our conclusions in section 
[VT| , where we also examine the implications for mf°' e . In Appendix [D] we point out simple 
extensions of either the MSSM or the GUT in which our (constrained) analysis does not 
apply. 

Below, we point out that (negative-energy) color-breaking local minima of the full scalar 
potential are generic. Approximate analytic constraints that are often used in the literature 
attempt to eliminate not only global but also local color-breaking minima, and are thus too 
strong. Furthermore, they are designed to eliminate only certain types of color-breaking 
minima, and are thus too weak. This is explicitly demonstrated in Appendix |B|, where we 
also compare analytic and numerical treatments of the problem. 

We follow Ref. || in calculating the couplings (and the unification point), and Ref. 



21] in treating the one- loop effective potential correction AV, including contributions from 



all sectors. The boundary conditions at M G for < m < 500 GeV, \A \ < 3m , and 
50 < Mi < 350 GeV, are picked at random. m^ ole values are picked at random but with a 

gaussian distribution defined by (^). For a given vn? t ole , tan (3(Mz) is picked at random within 
its allowed region [case (1)]. Values of m\ ole < 155 GeV are discarded so that tan/3 > 1.1 
and the divergent limit is not reached (see below). We then solve iteratively for (and 
Bq). [The relevant renormalization group equations (RGE's) have been given by various 
authors [^2j.] The sign of /io, which is a RG invariant, is picked at random as well. We 
explicitly verify that the potential is bounded from below; SU(2) x U(l) is properly broken; 
the lightest supersymmetric particle (LSP) is neutral (in practice, it is a neutralino); and 
that all the squared masses which correspond to physical scalars are positive before a point 
is accepted. We then find, in general, that we are in agreement with lower bounds coming 
from direct searches, (e.g., values of Mi < 100 GeV, which would imply a too light chargino, 
are already excluded in our case by a tachionic t-scalar.) We do not impose a lower bound 
on m h o, which is the subject of our investigation. Our treatment of color and/or charge 
breaking (CCB) minima is discussed in Appendix |B[ 
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II. THE WEAK-SCALE HIGGS SECTOR 



The Higgs part of the MSSM (weak-scale) scalar potential reads P3| , |TE 
V(H U H 2 ) 



(m 



-Bii{R x E 2 + h.c.) + 



Hl + u, 2 )\H 1 \ 2 + (m 2 H2 + u, 2 )\H 2 \ 

^MSSM 



-(\H 2 \ 2 - {H^Y + AV, 



(7) 



where m% l , m\ , and B (//) are the soft (supersymmetric) mass parameters renormalized 

down to the weak scale, m\ = Bjj, < 0, A MSSM = 92+ ^ 91 ; and we suppress 577(2) indices. 
The one-loop correction p4| , [l8| AV^ = /\v Aorie_/oop (which, in fact, is a threshold correction 
to the one-loop improved tree-level potential) can be absorbed to a good approximation in 
redefinitions of the tree-level parameters |H]j25[l . 

A broken SU(2) x U(l) (along with the constraint m 2 H +m 2 H2 + 2u 2 > 2\m\\ from vacuum 
stability) requires |23| , |T8H 



(m Hl + n ){m H2 + fi ) < |m 3 



2 1 2 



and the minimization conditions then give [23,18 



ni 



I'' 



Hi 



m 



Hi 



tan 2 13 1 



tan 2 (3 



-M 



(8) 



(9a) 



ran 



sin 2(3 



m 



(9b) 



For tan/3 — > 1 one has h t 3> /lb, 



and hence cannot have m 2 Hl 



m 



H 2 



(assuming universal 



initial conditions). Thus, —>■ oo in that limit [case (1)], and the SU{2) x U(l) breaking is 
driven by the B/j, term. In practice, we will stay away from the divergent limit in case (1) by 
taking tan/3 ^ 1.1 [21|. We do not expect significantly different results for 1 < tan/3 < 1.1; 
i.e., the divergence is either stabilized by model dependent finite-loop corrections or would 
exclude that region (see also below). For tan/3 — > oo [case (2)] one has B\i — > so that the 
symmetry breaking is driven by m 2 H2 < 0. 

The tan/3 — > 1 case corresponds to an approximate SU{2) L+R custodial symmetry of 
the vacuum p6|,^7l, which we further discuss in Appendix |A|. The symmetry is broken at 
the loop level so that one expects tan (3 slightly above unity, in agreement with our cut 
tan/3 ^1.1. As a result of the symmetry, the CP-even Higgs mass matrix becomes 



M 2 w 



/j 2 x 



-1 1 



(10) 



and it has a massless tree- level eigenvalue, 
of the tree-level formula [p3| , |18 l 



h° 



0. This is, of course, a well known result 



when taking /3 = |. The mass is then determined by the loop correction m h o m A h o oc h t m t 
and is further studied in sections and |V|. On the other hand, in case (2) (3 — > | and we have 
Mz (assuming m A o > M z ). When adding the loop correction^ m h o < V2M Z ~ 130 



m 



T 
h° 



GeV. The heavier CP-even Higgs boson mass eigenvalue equals approximately \/2|/i| in case 
(1) and rriHO ~ m A o in case (2). The loop corrections are less relevant here as typically 
m 2 H0 3> A^ (in particular, when \n\ — > oo). 

The custodial symmetry (or the large /i parameter) dictates in case (1) a degeneracy 
m A o »s mud ~ ~ v2|/jt| . [The tree-level corrections to that relation are of order 

(Mw, z/itia ) 2 -] That is, at a scale Q rs \/2|/x| ~ 2 TeV the heavy Higgs doublet H is 
decoupled, and the effective field theory below that scale has only one massless (m^ ~ 0) 
Coleman- Weinberg SM-like (u h o = v) Higgs doublet, h. Unlike in our case, in general 
M z /m A o is not guaranteed to be small, and the decoupling scale Q is much lower, typically 
below the scalar-quark thresholds. 

The Higgs sector in case (2) exhibits an approximate O4 x O4 symmetry (see Appendix 
|A|). This is a special case of the m| — » case discussed in Ref. |28|]. To obtain m| ~ 



(starting with universal initial conditions and hj, ~ h t ) one would require large Mi and 



|/xo|, which enable one to obtain the desired limit by adjusting large cancellations |29|-|3T 
Once expectation values are acquired, A and H + are massive pseudo-Goldstone bosons 
of the broken O4 symmetry. However, they are typically heavy (for Vdown = the O4 
symmetry is restored) but lighter than the scalar quarks in this scenario (whose mass rriq ~ 



ml + 6M? w 1 TeV). 



Large tan/3 [case (2)] solutions are discussed in Ref. p9| , |30| , p!3| , |32| , |3T| . They predict 
m,f 0le ~ 180 ± 15 GeV |pl,|T3| , heavy scalars, and contain potentially troublesome loop cor- 
rections, e.g., to m b , b — >• 57 [^3| , p2| ; and, in the absence of light scalars to counterbalance 
large m t contributions, also to Z — > bb (see Ref. |[33||). The former can eliminate the motiva- 
tion for any Yukawa coupling analysis. A different approach was recently presented in Ref. 
pl| where those large finite-loop corrections were used to construct models consistent with 
universality, a broken SU(2) x U(l), and h t = h b = h T . Finite superpartner (~ 1 TeV) loops 
are adjusted to diminish m& by the right amount (~ 1 GeV), and mb is strongly dependent 
on the values of the soft parameters, a possible but undesirable situation. Much lower values 
of m v t ole are then preferred. Diminishing m t may rectify the situation with Z — > bb, but the 
prediction for the b — > 57 branching ratio can be significantly modified by the same super- 
partner loops flT5| . In contrast, the tan/5 w 1 [case (1)] solution constrains only m v t ole ^ 140 
GeV {jr? t ole > 155 GeV when imposing the tan/5 > 1.1 cut), is in better agreement with (^j), 
and predicts a light Higgs boson {m h o < 110 GeV). Cases (1) and (2) are summarized and 
compared in Table |. We hereafter study the tan/5 w 1 [case (1)] solution only. We assume 
that finite-loop corrections to m b are at the most a few percent which is, in general, the case 
away from the large tan/3 limit, and thus do not alter our results. [The parameter space is 
sensitive to those corrections only in case (2).] 



'A/jO is enhanced in this case by the large mt and the heavy t-scalars. 
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III. THE tan/3 -► 1 SCENARIO 

Our primary motivation to study the tan /3 ~ 1 solution to the MSSM is its consistency 
with condition (fj). The large values of a s (M z ) implied by ([[]) [e.g., Eq. generate (too) 
large positive loop corrections to h b /h T , which need to be counterbalanced by large negative 
corrections oc hf, ti^. The latter are sufficient only if some of the the Yukawa couplings are 
large, i.e., near their unitarity upper bounds, which is the case in the allowed regions. [Also, 
H ag in (Q) becomes non- negligible in that limit.] The large values of the Yukawa couplings 
can be understood in terms of a quasi-fixed point in the flow of the respective RGE's |34 



In particular, in case (1) h t is near its quasi-fixed point. The two issues, i.e., the consistency 
of tan (3 ~ 1 with (fj) and the /i r RGE flow structure, are strongly related. 

The smallness of A^ (or alternatively, the large can be understood in terms of an 
approximate SU{2)l x SU{2)r symmetry in the Higgs potential, which, if it exists at some 
scale, will be only slightly broken at Mz, as is shown in Appendix [S]. The tan/5 ~ 1 



scenario was previously studied in Ref. ||35|| , where it was referred to as "highest classical 
degeneracy" : in the tan (3 = 1 limit condition @ becomes an equality. [In practice, we find 
that {m 2 Hl + /i 2 )(m^ 2 + /i 2 )/m| ^ 0.99 for A/3 < 1 and after the proper redefinitions.] It 



was subsequently studied in Ref. |36|] in the context of a global MSSM (i.e., no high-scale 



assumptions), and its consistency with limits on m^o was shown. The tan/3 « 1 solution 
suggested by the minimal SGUT relation = h T revived the interest in that corner of 
parameter space. More recently Ref. P^-B3, motivated by the quasi-fixed point prediction 



to m t , extended the discussion from the Higgs sector to the full parameter space. 

In Figs. |]a an d b we show the prediction for \i as a function of m^ ole and of tan/3, 
respectively. Typically ~ 1 TeV, depending on A^ and on the soft parameter scale. It is 
the large parameter that dictates the characteristics of the scenario (and not only for the 
Higgs sector). Before further elaborating on the calculation of m h o, let us briefly summarize 
some of the features that appear in our numerical studies, and which are of relevance for the 
discussion below. Though we agree on the characteristics with other authors, our emphasis 
and interpretation are different. 

1. The Higgsino (with mass ~ fi) is decoupled from the (much lighter) gauginos. The 
bino and the wino are the lightest neutralino and chargino, respectively. (The former 
is the LSP which has a bino fraction near unity.) The heavy Higgsino decouples from 
the Z — > bb vertex, and thus, smaller values of m? t ole are favored^. However, rrP° le 
and tan/3 are correlated, and as vnP t ole grows is diminished (see Fig. |2|). Thus, in 
principle, Higgsino - t-scalar loops can still counterbalance large m t contributions to 
Z — ► bb (see, for example, Ref. |33]|) if the t-scalar is light enough. [A light t-scalar is, 



however, less likely if \n\ is small (see below) and from the p-parameter.] 
2. The t-scalar mass-squared matrix is (we take all parameters to be real) 



6 The prediction for the Z — > bb branching ratio decreases with mf oe , and its measured value is 
~ 1.5 s.d. higher than the SM prediction with mf°' e ~ 140 GeV. The discrepancy grows with m^ ole 



s 



m| ™>t(At + fJ>/ tan/3) 



m t (A t + fi/ tan/3) m 



2 



(12) 



and similarly for the superpartners of the other fermions. The diagonal elements 
correspond to the left and right-handed t-scalar squared masses, m~ L , which consist 

of soft, F, and D terms. A t is the trilinear soft parameter preceeding the /it-f^z^i? 
term in the scalar potential, with A t (Mc) = Aq. The prediction for the diagonal terms 
has no unique charachteristics in our case, and the left and right-handed masses are 
in the ~ 100 GeV —1 TeV range. However, the mixing term is typically large (unless 
the fermion mass is ~ as is the case for the light families). In the limit |/i| — > oo 
there is a 9t = \ left-right mixing so that one of the eigenstates would have a negative 
mass squared, and the limit cannot be fully realized. In practice, the soft terms and 
mt\\i\ are often of the same order of magnitude, with a nearly degenerate mass-squared 
matrix. Thus, the requirement that (§) defines a physical minimum with no negative 
squared masses, corresponding to physical scalars, constrains /i. Also, our previous 
comment regarding a light t-scalar and a small /i-parameter is clear by observation of 
(|12"D. It is interesting to note that by fine adjustments of the soft parameters [i.e., of 
the degeneracy in flT2"D] one may realize a scenario with a light (< 45 GeV) t-scalar 
which is nearly decoupled from the Z (e.g., sin 6^ « 0.8). (See, for example, Ref. |[41|| .) 
However, such scenarios are often associated with false vacua (see below and Appendix 
©■ 

3. The absence of a negative t-scalar squared mass determines a (mo and Ap dependent) 
lower bound on Mi. In particular, "no-scale" (i.e., m = A Q = 0) SGUT models are 

strongly constrained^ [|. Also, requiring that the light eigenstate of the r-scalar mass- 
squared matrix (which carries charge) is not the LSP determines an upper bound on 
Mi in the no-scale limit (see also Kane et al. |2I|). These point at a correlation between 



the bounds on the soft parameters. A different class of correlations which take effect 
in the h t — > h{ lxed limit is discussed by Carena et al. [37]. Strong correlations (e.g. 



the mP t ole — tan /3 correlation^, the h t h{ txed ones, and those due to scalar mixings 
discussed here) do not allow elimination of the large- solutions (that correspond to 
large left-right mixing enhancement to m h o - see below) on the basis of fine-tuning 
arguments. The naive notion of fine-tuning becomes obscure and misleading in the 
presence of strong correlations. (See also Ref. []3T|].) Rather, strong correlations make 
the scenario quite predictive. 

4. We can rewrite the mixing term in ( |T2"D as htH^T where = v uv . f is strongly 
constrained by requiring that the full scalar potential does not have a color and/or 



7 The strict no-scale assumption Bq = uiq = Aq = is not consistent here. See Appendix 



No-scale models with hf, = h T were recently discussed in Ref. [42] 



9 Note that even if we do not require hb = h T , mf ole and the lower bound on tan j3 are correlated 
due to unitarity considerations. 



9 



charge breaking global (GCCB) minimum, which is a different issue than whether the 
extremum defined by @ is indeed a minimum. The issue of false vacua is of particular 
interest in the tan/3 — > 1 scenario as V(Hi, H 2 ) — > 0~ in that limit. Below, we will 
consider that subject in some detail. In particular, see Appendix [B]. We find that the 
absence of negative energy CCB minima implies the constraint, first described in Ref. 

m, 

(\A t \ + s\ t i\) 2 <2(ml+m 2 iR ), (13) 

where s = A t fi/\A t fi\, (\A t \ + s\fi\) 2 = f 2 (for tan/3 = 1), and m 2 ^ consist of only 
the soft terms. (Note the importance of the trilinear parameter and of the relative 
sign.) Our numerical studies imply that fll3|) is (to a good approximation) a sufficient 
condition (in that region of parameter space) but is not necessary (i.e., it is too re- 
strictive). For example, more than 80% of the points that are inconsistent with that 
constraint for m v t ole > 165 GeV correspond to CCB minima which are only local and 
are therefore safe. In Fig. || only filled squares correspond to (probably unacceptable) 
GCBB minima. Also observe in Fig. |^ that fi < is preferred by CCB constraints. 
This is because typically A t ^> Aq for Aq < so that the contribution from a large 
and positive \i cannot be canceled in ( |T3|) (for our range of A ). 



IV. THE LOOP-INDUCED MASS 

Let us now proceed to discuss the mass of the SM-like higgs boson. The upper bound 



on rriho is given, e.g., in Ref. [P 
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and where m~ are the eigenvalues of the t-scalar mass-squared matrix, Q t is the mixing 
angle, and we have neglected other loop contributions. (They are included in our numerical 
calculations below.) The tree-level mass squared, m^ 2 , and the loop correction, A 2 ^, are 
bounded by the first and second terms on the r.h.s. of Eq. ([14]) , respectively. In the absence 
of mixings A$ t = 0. For tan/3 — > 1 one obtains m^ — > 0, and thus m^o pa A^o. 

The upper bound on A^o scales as m\ and depends on the overall scale of the t-scalar 
mass, the amount of left-right mixing, and the separation \m~ — mjj. The large mixing and 
separation for — > oo enhance A h o significantly (typically by 10 — 20 GeV) compared to 
the first term alone (which was all that was included in Ref. fi~2||). Condition (§) (and h t 
unitarity) correlate m\ ole and tan/3 in case (1) (see Ref. P,[16|), and A^ increases from ~ 
for mP t ole ~ 140 GeV to ~ 1 for m v t ole ~ 190 GeV. By increasing tan /3 one gradually departs 
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from the — > oo limit, and the enhancement to A^o from left-right mixing decreases. Thus, 
there is an interplay between Ag t and the overall factor of m\. (Of course, m^ grows with 
tan/? and hence with m pole .) 

As was pointed out in the previous section (see also Appendix |B|), \i and the amount 
of mixing T are subject to color and/or charge breaking (CCB) constraints. In this section 
points which are consistent with (|T3|) , i.e., which have no negative-energy CCB minima, are 
marked by a circle - O. Points which are inconsistent with (O) and therefore have either a 
local or global negative-energy CCB minimum are marked by a "diamond"- O. The filled 
diamonds correspond to a global color breaking minimum (as determined using numerical 
procedures), which is (most likely) unacceptable. The open diamonds represent CCB minima 
which are only local (i.e., above the standard model minimum) and are therefore safe. 
Constraint ( fL3|) is always evaluated at the t-scalar scale, independent of the scale chosen for 
minimization. The same is true for the numerical search of global CCB (GCCB) minima. 

Calculating m h o is reduced (in our case) to a great extent to calculating A^o, i.e., the 
deviation from what would be the case if the matter were supersymmetric at the weak scale 
[e.g., mf = mi 2 = m t in (O ) and flTa])]. Two standard ways to perform the calculation are 
the effective potential method (EPM) (e.g., see Ref. [[I5|,[l4],f|6|,[3(| ) and the renormalization 
group method (RGM) (e.g., see Ref. |^ , |48|]36|| ) . The two methods correspond to a "run and 
diagonalize" and "diagonalize and run" algorithm, respectively. The algorithms and their 
implementations are discussed and compared in greater detail in Appendix |Cj. Here we would 
like to demonstrate that the EP and RG methods do agree (within their domains of validity), 
i.e., the results are not sensitive to the way one treats the threshold corrections. In particular, 
both methods exhibit a large (10 — 20%) two-loop ambiguity from the presence of some two- 
loop residual terms in the calculations (next-to-leading and leading two-loop logarithms, as 
well as finite terms). We explicitly identify and discuss such terms in Appendix |CJ The 
residual logarithmic two-loop terms can be described in terms of a scale ambiguity (i.e., 
the scale dependence of the one- loop calculation is of two-loop order). Indeed, comparing 
the two methods and studying that ambiguity allow us to estimate the two-loop correction. 
However, the uncertainty remains. 

In the EPM the ambiguity is explicitly related to a scale dependence, as one is required to 
specify a subtraction scale Q (i.e., halt running at Q and diagonalize mass matrices). Indeed, 
including AV one ~ loop (Q) corrects for one-loop leading logarithms, and thus guarantees a 
scale independence of the calculation, but only to that order (and up to wave-function 
renormalization). It is suggestive to take Q = Mz, i.e., the scale at which the physical 
inputs are given. Another attractive choice is a scale at which AV two ~ loop (Q) is minimized. 
We will not attempt to find such a scale rigorously, but estimate Q ~ 600 GeV. (One indeed 
expects such a Q to be in the vicinity of some average scalar-quark scale, as for /\y one - loo v 
in Ref. We will use the comparison between the two methods to justify our choice. 

(Of course, there is no one choice that is suitable for all sets of parameters.) 

In the RGM, which sums the leading logarithms to all orders if one integrates numerically, 
the ambiguity is more apparent in terms of orders in perturbation theory (i.e., the order 
at which the summation is truncated in approximate iterative solutions). However, there is 
also an explicit scale dependence from the choice of decoupling scales for the heavy fields 
(i.e., decouple heavy eigenstates, calculate RGE's, and run to the next decoupling scale). 
We will decouple the heavy Higgs doublet H at m^o « 2.1 TeV and other fields at lower 
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thresholds^. When decoupling H we include appropriate matching conditions ||8 . 

The EPM is more straightforward to implement (and is therefore more common in the 
literature) while the RGM requires a careful consideration of decoupling scales and matching 
conditions. 

The scatter plot Fig. |3| displays the results of an EPM calculation with Q = Mz- (When 
increasing Q the soft mass parameters and m t decrease so that the EPM m^o is diminished.) 
The scatter plot Fig. ^ displays the results of a RGM calculation (performed numerically). 
Comparing them, one observes a 10 — 20 GeV difference between the respective Higgs boson 
masses (for a given point in the parameter space). However, this is not due to the different 
choices of an algorithm, but to the different residual two-loop dependences in the two calcu- 
lations. To illustrate the maximum ambiguity we used the t-quark pole (running) mass in 
the calculations presented in Fig. |3] (Fig. [|). (The ~ 5% difference in the t-quark masses is 
formally of two- loop order.) In both plots one can observe a separation into a heavier and 
a lighter Higgs boson branches. This is due to the enhancement to m h o in some cases from 
t L —tu mixing, and is further discussed in the following section. 

To further compare the two methods, it is useful to consider the following exercise. 
The RGM calculation could be modified so that leading logarithms are not summed to all 
orders, e.g., one could use the one- loop leading- logarithm approximation formula of Ref. 
17]. Alternatively, we could minimize the unsummed \a{mi/Q) leading (and/or next-to- 



leading) logarithms in the EPM by choosing Q m m^. In either case the residual terms in 
the different calculations are qualitatively brought on the same footing. We carried out that 
exercise numerically, and Fig. |5] (Fig. |6|) compare the EPM with Q = Mz {Q = 600 GeV) 
calculation to the RGM leading logarithm approximation (numerical calculation). Once the 
residual two- loop dependence of the calculation is properly adjusted, the two algorithms 
agree well, in particular for Q « 600 GeV. [All calculations presented in Fig. |5| (Fig. []) use 
the i-quark pole (running) mass.] 

This all points at the most important source of uncertainty: residual two-loop (and 
wave-function) scale dependences and a ~ 5% a s dependence of m t . Two-loop terms fuel 
the ambiguity, which therefore cannot be removed in a satisfactory fashion within a one-loop 
calculation. (Even though we could identify some two-loop terms, the consistency of their 
removal is doubtful). Expanding the calculation to higher orders is not straightforward. We 
discuss two-loop calculations in Appendix |CJ 

Figs. ^| and [| define the theoretical ambiguity in determining the loop induced m^o in 
this scenario. While the former sets the upper bound (for a given mf oZe ), the latter is 
a more accurate leading-logarithm calculation (and uses properly the running t-quark 
mass). However, within a one- loop calculation, and in particular when next-to- leading 
logarithms and some finite terms are comparable to higher-loop leading logarithms (which 
is the case here - see Appendix 0), neither method provides a more consistent calculation 
then the other. As we argue in Appendix |Cj, the former is also an upper bound on two- 
loop calculations, which most probably diminish m^o. We do not think that any stronger 
statements are justified at present. We continue to discuss the prediction and its upper 



10 We checked that the details of the choices for those intermediate scales introduce a ~ 1% 
uncertainty. 
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bound in the next section. 

A different issue is that of the constraints on the t-scalar left-right mixing enhancement to 
A/jO. In the EPM [e.g., Eq. (0)] the enhancement is a straightforward result of diagonalizing 
the CP-even mass matrix. In the RGM one has to define the decoupling scale for the t- 
scalars. If mf 2 then one would decouple only t\ while £2 loops could still contribute to 

A/jO. (In practice, the effect is accounted for in the RGM by a set of appropriate matching 
conditions.) However, the amount of mixing T (and thus the enhancement) is strongly 
constrained by the physical vacuum. 

Requiring that the SM extremum is a minimum (i.e., no negative squared masses cor- 
responding to physical degrees of freedom), and furthermore, is the global minimum (i.e., 
there is no deeper minimum that does not conserve color and/or charge) strongly constrains 
T. Constraints imposed by the latter may be evadedP] if the non-SM global minimum is 
separated from the SM local minimum by a tunneling time greater than the age of the 
universe ||50|| . Constraint (|I3|) (naively) attempts to eliminate all, i.e., local and global, 
negative-energy CCB minima in that region of parameter space and is too strong. The 
points marked in the scatter plots by only a an open "diamond" survive a more careful 
analysis, i.e., they correspond to only local minimaQ. We find, using numerical methods, 
that for m p t ole ^ 165 GeV (165 GeV < mf e ) roughly 80% (15%) of the points that fail to 
satisfy constraint (|H|) correspond to (probably unacceptable) GCCB minima of the (one- 
loop improved) tree-level scalar potential. The enhancement is then partially washed away 
by requiring consistency with the physical vacuum. The possibilty of coexistence of the two 
vacua (with a tunneling time greater than the age of the universe) is, at the present stage of 
our calculation, a remaining ambiguity. We further discuss those issues in the next section 
and, in greater detail, in Appendix |B[ 



V. THE PREDICTION FOR THE HIGGS BOSON MASS AND ITS UPPER 

BOUND 

In this section we present and discuss our numerical results. In section [IV] we discussed 
various methods to calculate the mass (to one- loop). The effective potential method (EPM) 
with a subtraction scale Q = Mz (e.g., see Fig. |3|) puts the weaker constraints on m^o and is 
the conservative choice for the upper bound. In Fig. [7] we present monte-carlo distributions 
of mho using that method for 155 < m? t ole < 165, 165 < mF t o[e < 175, and 175 < vrt? t ole < 185 
GeV. The upper bound can be read off the figures, i.e., m h o < 101, 106, 114 GeV for m? t ole < 
165, 175, 185 GeV, respectively. (Note that these bounds are deduced from our monte-carlo 
calculations, and are not rigorous.) Nevertheless, we would like to stress that typically m h o 
lies below its upper bound. Indeed, one cannot rule out the possibility (if mF° le is heavy 



This claim, however, does not take into account finite-temperature effects. 

12 In fact, all points are found to have at least a negative-energy CCB local minimum of type 
M 2 = - see Appendix [b|. However, we find that in our region of the parameter space (|i~3| ) is a 
useful tool to identify minima that compete with the standard model minimum. 
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enough) that is just outside the range that LEPII may cover. However, this is not likely 
on both grounds, the mF t ole distribution [Eq. @] and the m h o (theoretical) distribution (as 
reflected in Fig. |7|). 

The role of the large left-right mixing in the i-scalar sector can be seen in Figs. [5]a and 
|7|b, which exhibit two peaks in the mass distribution, a "logarithmic" and a "mixing" peak. 
The former is due to the diagonal entries in the t-scalar mass matrix [Eq. (|12D1 while the 
latter, which is at larger mass, is due to the mixing enhancement [Ag t in Eq. flip])]. The 
interplay between the mixing and the m\ factor in the loop correction A/jO [e.g., in Eq. ([14]) 1, 
and between A h o and the tree-level mass m[o (that grows with tan/3 and, therefore, with 
mf°' e ) can be seen by comparing Figs. [7|a-b with Fig. 0c. The mixing enhancement plays an 
important role for smaller values of mf oZe (A/? — > 0). In fact, that interplay generates for 
a fixed set of soft parameters a local minimum in the m^o distribution (not shown in the 
figures) near m? t ole ^165 GeV. 

The left-right mixing enhancement is strongly constrained when requiring that the min- 
imum of the full scalar potential is the physical one, i.e., no fields other than H\ = Vdown 
and H® = v uv have non- vanishing expectation values (see section [I I] and Appendix ]B|). 
Constraint ( |13|) eliminates potentially dangerous points in the parameter space (marked by 
a "diamond" in Figs. and |j), some of which, however, correspond to only a local color- 
breaking minimum (marked by an open "diamond"). Those which correspond to global 
minima are marked by a filled "diamond" . Fig. |3| is the same as Fig. |7| except that points 
which correspond to global color and/or charge breaking (GCCB) minima are omitted (the 
shaded areas have no negative-energy CCB minimum, global or local.) The mixing peak 
has diminished significantly but has not disappeared. Its presence is even stronger for larger 
values of mf le (see appendix ||). (Some caution is required: the cosmological history of the 
universe may accommodate such a global minimum.) 

A different issue is that of the large two-loop ambiguity in the calculation pointed out 
in section [TVJ. For example, in Fig. |S| we show the same distributions as in Fig. [7], but in 
which the running (rather than the pole) t-quark mass was used to calculate m/jo. (The 
running mass is ~ 5% smaller.) Also, the distribution corresponding to the renormalization 
group method (RGM) calculation (e.g., see Fig. £|), which we showed^ in Fig. || to roughly 
agree with the EPM calculation (with Q « m ( -), is shown in Fig. |10] for 155 < m p t ole < 165 
and 165 < mf°' e < 175, respectively. Again, the t-quark running mass is used. Also shown 
are the distributions when GCCB points are omitted (shaded areas) to be compared with 
the total areas (shaded and unshaded) in Fig. [8] Note the concentration at lighter masses 
compared to previous figures, even though the upper limits are only slightly changed. It 
should be stressed that had we not included the proper RGM matching functions pS| (which 
would roughly correspond to setting the left-right mixings to zero by hand) then the RGM 
distribution would sharply peak at the lower end of the m^o range. Thus, the importance 
of the matching functions is obvious. 

All, Figs. - [HI are consistent within the theoretical ambiguity. It is suggestive that a 



The RGM calculation sums the leading logarithms to all orders in perturbation theory, while 
in the EPM calculation with Q ~ higher order terms are, in principle, minimized, and do not 
need to be summed. 
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consistent two-loop calculation will reduce the upper bounds given by the EPM calculation 
(we choose conservatively Fig. |?]), but probably not to lower values than suggested by the 
RGM (Fig. [10]). Further study is required before any stronger conclusions can be drawn. 
Even so, for mf°' e < 175 GeV m h o [in case (1)] is within the reach of LEPII, and possibly 
of the Tevatron (see Ref. |5TJ). If m p t ole is lighter, as suggested by electroweak precision 
data, our histograms imply that m/jo may still be relevant for Higgs searches at LEPI. (In 
particular, if two- loop corrections reduce m h o as suggested by Fig. |TTj.) 



VI. CONCLUSIONS 

To summarize, we worked in a constrained framework assuming (up to high and low-scale 
threshold corrections) a>i = a 2 = = ac, h b = h T , and universal initial conditions. Also, 
we explicitly included the correlation between the t-quark mass and the weak angle from 
electroweak precision data. Within that framework, only tan (3 ~ 1 and large tan /3 regions 
are allowed. In the large tan/3 region (which is allowed only for a large mf oZe , unless nib 
is allowed to get large finite corrections from superpartner loops) m^o is bounded by the 
triviality bound, and is probably not interesting for Higgs searches at LEPII. 

In the tan /3 ~ 1 region the tree- level rrv[ nearly vanishes, but grows with m\ ole . The loop 
correction to the mass is enhanced by large t^ — tn mixing terms. The latter are diminished 
with an increasing m pole but, as m? t ole grows, the overall factor of m\ can compensate for 
that effect. The enhancement due to the large mixing is further constrained by requiring 
the physical vacuum to correspond to the global minimum of the full scalar potential. We 
find that for parameters corresponding to large mixing there is usually a color breaking 
minimum, but in many cases this is only a (presumably harmless) local minimum. We also 
pointed out that a two-loop ambiguity in the calculation of the Higgs boson mass is enhanced 
in our case and can be as large as 20 GeV. All these effects lead to the conservative bound 
m h o < 100 (110) GeV for m? t ole < 160 (175) GeV. (We briefly discuss two-loop calculations 
in Appendix y, where we conclude that such a calculation will most probably strengthen 
the upper bound.) The Higgs boson mass in the class of models characterized by tan/3 ~ 1 
is then almost certainly within the reach of LEPII (and possibly within the reach of the 
Tevatron ^TJ). Furthermore, the Higgs boson is typically lighter than the upper bound, and 
may be light enough to still be discovered at LEPI. 

Regarding the t-quark mass, one requires for SU(2) x U(l) breaking (in this scenario) 
m poie > ~l4q Q e y ; anc } f ur ther requiring tan/3 7^ 1.1 (which is the case if the divergent limit 
is to be avoided, i.e., if the divergence is not stabilized by model-dependent finite terms) we 
obtain m? t ole > 155 GeV. Eq. (|5|) and, in particular, the structure of the radiative corrections 
to the Z — > bb vertex, however, suggest that rrP t ole is not much larger. (The latter can be 
evaded only if both the t-scalar and the Higgsino are light, which is not very likely in our 
case.) Also, we were able to understand the smallness of the parameter space allowed by 
hi, = h T in terms of a custodial symmetry. The symmetry (theoretical) argument suggests 
m poie < Q e y < 1). In Appendix [D]we briefly examine extended Higgs sectors, and 
point out that the Higgs mass is a possible probe of such extensions if mF t ole (and/or tan/3) 
is known. 

In Appendix |B| we point out that the full (one- loop improved) scalar potential generically 
possesses color and/or charge breaking negative-energy local minima. Analytic constraints 
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that are often found in the literature attempt to eliminate all such local minima (i.e., are 
too strong) but by construction apply to only certain types of minima while ignoring others 
(i.e., are too weak). We give a general classification of the local minima and are able to 
further identify [in case (1)] those that can potentially compete with the standard model 
minimum. 
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APPENDIX A: THE HIGGS SECTOR CUSTODIAL SYMMETRIES 

In this appendix we examine some symmetry aspects of the SGUT solutions, and in 
particular, of case (1). Some features of the latter, as well as the distinction between 
the two cases, can be better understood in terms of the custodial SU{2)l+r approximate 




symmetry of the SM lagrangian [pqj ; i.e., turning off hypercharge and flavor mixings, and 
if h t = h b = h, then one can rewrite the t and b Yukawa terms as an SU{2) L x SU{2) R 
invariant, i.e., 

(where in the SM H 2 = ir 2 Hf) and SU(2) L x SU(2) R -> SU(2) L+R for v ^ 0. However, 
h t ^ hb and the different hypercharges of t c L and b c L explicitly break the left-right symmetry, 
and therefore the residual custodial symmetry. 

In the MSSM, on the other hand, Hi is distinct from H 2 and if v down ^ v up SU(2)i x 
SU(2) R — ► U(1)t 3L +t 3R p7j - In particular, the SU(2)l+r symmetry is preserved if 8 = | 
but is maximally broken if 8 — ^. Ignoring hypercharge, we can then define for illustration 
two quantities, 

hf — hb 

Bm = v down , (A2a) 

m t - m b 



„ tan 8 — 1 , . _ . 

B H = m b , A2b 

m t - m b 

that measure the SU(2)l+r symmetry breaking in the matter and Higgs sectors, respectively. 
If we map the tan/3 — mF t ole plane onto the Bm — Bh plane, we find that cases (1) and (2) 
approximately correspond to (1,0) and (0,1), respectively. Thus, while the constrained 
5*0(10) solution realizes m t ^ nib as an intrinsic property of the vacuum (i.e., v up ^> v down ) 
and the grand-unified symmetry protects Bm ~ |29j, the tan/? w 1 case realizes the 
SU(2)l+r (explicit) breaking as a scaling phenomenon (i.e., Bm <x h{ lxed ) and the vacuum 
respects the symmetry (up to loop corrections). 
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The symmetry that is spontaneously broken in case (2) is a O4 x O4 symmetry. Haber and 



Pomarol observed [28] that for m 3 — > [which is the situation in case (2)] there is no mixing 
between H 1 and H 2 and the Higgs sector respects 4 x 4 (up to gauge-coupling corrections). 
The symmetry is broken to 3 x 3 for v dmun 7^ v up 7^ and the six Goldstone bosons are the 
three SM Goldstone bosons, A , and H ± . The symmetry is explicitly broken for g 2 7^ (so 
rriH+ = My/) and is not exact even when neglecting gauge couplings (i.e., m 3 7^ 0). Thus, A 
and H ± are massive pseudo- Goldstone bosons, i.e., m 2 H+ — ~ m 2 A0 = C x m|. However, 
C = — sin 2 2 ^ and can be large, which is a manifestation of the fact that O4 x 4 -> O4 x O3 
for ^own = 0. (The limit m 3 — > corresponds also to a £7(1) Peccei-Quinn symmetry JIB].) 

Our main interest in this work is, however, case (1), where from fl9b| ) the Higgs-sector 
SU (2) custodial symmetry constrains 

2,22,2 2 / A Q\ 

m Hl + fx = m H2 +fi = -m 3 . (A3) 

This is realized by taking the limit |/i| — > 00, provided \B\ « which in our formalism 
comes out as a prediction. (Indeed, many SGUT models can easily realize a large B , 
for example, see Ref. [0].) Relation (|A3|) results in a (tree-level) massless eigenvalue in 
the neutral CP-even mass matrix. (Note that the SM Goldstone bosons are also massless 
Goldstone bosons of the broken SU(2) x SU(2).) 

The smallness of is protected by the custodial symmetry, i.e., 7^ is explained 
by loop corrections. The broken symmetry in the Yukawa sector induces at the loop level a 
small breaking in the Higgs sector (e.g., m 2 H 7^ rn 2 H2 ). The magnitude of the effect can be 
estimated by integrating c?tan/3/ tan/3, i.e., 

A p (M z ) = Ap(M G ) -JL£^(h 2 t - h 2 b )d\n^- » Ap(M G ) + 0.6, (A4) 

where we neglected h T , as well as threshold corrections, took h t w 1, hj, « 0, and assumed 
A^ < 1. Thus, even though the parameter space in case (1) is severely restricted, we can 
understand the restriction in terms of a symmetry. 



APPENDIX B: CONSTRAINTS ON COLOR AND CHARGE BREAKING 

MINIMA 

Let us consider the super-potential 

W = h t t L t R H° 2 + pH°H%, (Bl) 

where we have performed an 577(2) rotation so that H 2 has no vacuum expectation value. 
(We do not distinguish our notation for a superfield from that of its scalar component.) It 
is easy to convince oneself that only those terms in the super-potential are relevant when 
searching for the global minimum of the scalar potential in the large h t limit ||53|| . Using 
standard techniques we arrive at the corresponding scalar potential: 

h 2 V = M 2 — T + A, (B2) 

where the bilinear, trilinear, and quartic terms are 
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M 2 = m\Hf + m\Hf + 2m 2 # 1 °# 2 ° + m\? L + m\ R P Ri 



(B3) 



2(\A t \H 2 +s\fi\H°)t L t R 



and 



A = (t 2 L + i 2 R )Hf + P L P R + ^(- x U D - terms"), 



1 ,7T 



(B4) 



(B5) 



respectively. All fields were scaled 



4>/h t and are taken to be real and positive (our 



phase choice for the fields, which fixed m\ < and T > 0, i.e., maximized the negative 
contributions to V) and all parameters are real, m 2 2 = rn 2 Hi „ + /-t 2 , s = fiA t /\fiA t \, and the 



expression for the 11 D— terms" is |53 
12 



-Oil 



ET Ef t\ 2t 



+ 



2 f 6 3 



R 



+ a 2 



nf - Hi" + 1 



2 , J2 



p _ 72 



(B6) 



We can parametrize the four fields in terms of an overall scale X and three angles < 
a, P, 7 < |: = X sin a cos /3, E® = X sin a sin /3, tj? = X cos a cos 7, i/, = X cos a sin 7, 
and redefine 



h 2 V(X) = M 2 (a, p, 7 )X 2 - r(a, /?, 7 )X 3 + A(a, /3, 7 )X 4 . 



(B7) 



Then, for fixed angles, V(X) will have a minimum for X 7^ provided the condition 
32M 2 A < 9r 2 is satisfied. In that case, 



Y 



3_r 

8A 



and 



h 2 V ■ 

1 v mm 



'-XL 



32M 2 A^ 

9r 2 



> 0, 



(B8) 



M' 



(B9) 



We only allow parameters for which a standard model (SM) minimum exists. The SM 

/ "down is the 



up 



minimum corresponds to a = | and /3 = (3 (7 is irrelevant and tan /?° 
angle used to fix /1, m 3 , as well as the Yukawa couplings). It is easy to convince oneself that 
in that limit the 4x4 second-derivative matrix is 2 x 2 block diagonal (otherwise baryon 
number is violated). Thus, it is sufficient to confirm that the four physical eigenvalues are 
positive to ensure that it is a minimum. (This is done in our numerical calculations.) If these 
conditions are satisfied then the SM is at least a local (negative-energy) minimum, and one 
hasr SM - n »' 2 

'H v min — M SM/ 4il 



0, M| M < 0, and 



8D and 



reduce to the usual results X™ n 



-M 2 M /2A 



SM- 



rat 



Let us now consider the possibility of additional CCB minima with cos a 7^ 0. Since 
> for the class of models we are considering this requires f ^ 0, which we assume 



hereafter. From ( |B7| ) - ( |B9| ) it is easy to classify the possible color and/or charge breaking 
(CCB) minima for definite a, 0, 7. One finds that for T 2 < 32AM 2 /9 there is no CCB 
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minimum, while for 32AM 2 /9 < T 2 < 4AM 2 the CCB minimum exists but has V^f n B > > 
V™ n , which is presumably safe. For 4AM 2 < T 2 (including the more rare case M 2 < 0, that 
must fall in this category) there is a negative-value CCB minimum, which may however be 
either local (presumably safe), i.e., V£f£ > Knin; or global (probably unacceptable), i.e., 
Vmfn < VmL- Here, a sufficient (but not necessary) condition for an acceptable model is 

T 2 < 4AM 2 , (BIO) 

which generalizes that of Ref . |)3] . In principle, the above discussion holds for any number 
of fields (i.e., any number of angles), only the explicit expressions for M 2 , T, and A are 
more complicated. If constraint ( |B10[ ) holds for every choice of a, (3, and 7 (cos a 7^ 0) 
then there is no negative- valued color and/or charge breaking minimum, global (GCCB) or 
local. In the special case M 2 < the constraint cannot be satisfied and there will be a 



negative-energy local CCB minimum, an observation made previously by Gunion et al. [53 



If ( |B10|) does not hold, further investigation is needed to determine whether the minimum 



is global or local. It is straightforward to show that ( |B10| ) is satisfied for sin (3 = 0, i.e., there 



is no negative-valued CCB minimum for ifg = f°r the class of models we are considering. 
We will therefore restrict our attention to the case 7^ 0, in which case it is convenient 
to reparametrize X = an d rescale all fields by an additional factor of l/H®- This 
simplifies the discussion and enables us to examine and generalize previous work on the 
subject @|||5|,|3|. 

By inspection ( |B10| ) cannot hold for all values of the fields and there is always (and not 
just for M 2 < 0, which is hard to achieve in dangerous directions) a local CCB minimum, 
e.g., for M 2 = (A is positive definite). For example, one can choose 

Hi = -^f, (Blla) 
mf 



t L 



\2 



l(m \-mi], (Bllb) 



m 



~ ~ 7Tt~ 

i R = h-^, (Bile) 



mi 



which is always possible because of (J8|). We confirmed this in our numerical studies, i.e., 
there is a CCB (negative-valued) local minimum even if V SM = V(v down , u up ) is the global 
minimum. (We do not consider additive constants.) Although such local minima may be of 
cosmo logical interest, the relevant question in our case is whether V SM is the global minimum 
of the full (or more precisely, of the four- field) scalar potential. Similarly, we assume that 
positive- valued CCB minima would not be populated and are therefore harmless. The answer 
has to be given by numerical mapping of all minima, but it is still useful to review how one 
could derive analytic constraints from (|B10|) , which are typically relevant only for specific 
regions of the parameter space. 

For example, if we fix Hf = 0,t L = Ir = 1, then ( |B10|) gives the well known result [Q 

A\ < 3(m? L + ml + m\ 2 + ^ 2 ). (B12) 
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The equal field direction was chosen so the positive D— term contribution oc 1/h 2 vanishes. 
This is not a relevant requirement for h t Pd 1 (i.e., there may exist a deeper minimum with 
non vanishing D— terms). It is obvious that ( |B12| ) is not relevant when — > oo. We find 



that a more useful constraint is [taking Hf ~l,ti~t^ = t^Cl and using (|B10|) 1 



<2(mf L +ml R ] 



(B13) 



Note that if the order t 2 corrections to (P13|) are not negligible then V is less negative, which 
motivated our choice t < 1, i.e., the more dangerous direction. That constraint was also 
used by Drees et al. [[f3[] to describe their numerical results (i.e., if we take rri\ = m 2 in their 
formula). We confirmed in our numerical studies that (p3 1 3|) and not ( |B12|) is relevant in 
our case. [We found nearly all points to be consistent with (|B12|) for our ranges of the soft 
parameters.] 

Finally, we perform numerical minimization of the potential (|B~2|) and map all minima 



using monte-carlo routines. We find that (|B13j ) is to a good approximation sufficient to 
avoid a GCCB minimum, but is not necessary. For example, only about 80% (15%) of 
the points which are inconsistent with ( P13|) correspond to GCCB for m? t ole < 165 GeV 
(m^ ole ^ 165 GeV). We compare the points consistent with ( |B13j ) and with the numerical 
search in Fig. |[ The greater relevance of ( |B13|) for m v t ole < 165 GeV can be understood if we 
recall that the lower bound on | cos2/?| increases with mf oie , and thus V SM becomes a deeper 
minimum. (The lower bound is related to h t unitarity, and that statement is independent 
of our hf, = h T assumption.) On the other hand, and the negative terms in V CCB [e.g., 



That interplay shifts the CCB minimum from a 



Eq. ([B2|) 1 diminish with increasing m v ° le 
global to a local minimum. 

Let us demonstrate that claim in a simple-minded example assuming h t m h{ 
We define sin 2 /3 = 2m 2 /h 2 v 2 = x, 0.5 < x < 1 [u — 2M w /g 2 ]. We can rewrite 



ixed 



1-1.1. 



h 2 V SM 



-M 2 z h 2 t v 2 [1 - 4x(l - x)\ . 



(B14) 



|/i2y SM | sharply increases with m\ ole (i.e., with x). On the other hand, taking m 2 H2 w M| ~ 
in (Ral) we have 



4 



1 — X 
2x-l 



(B15) 



and fi 2 decreases from +oo to in the above range of x. Now assume that constraint ( |B13| ) 
is not satisfied, i.e., 



2t aCCB 



1 

— J 

2 



■ — Kfi [K/J, 



m 



tL 



mi) t\ 



(B16) 



where t 4 <C t 2 and k is a positive number in the range [|, |] (X min = K\fi\). We also assumed 
\A t \ <C and neglected t A terms. (M 2 < would increase k.) To further simplify, consider 



the case m = 0, so that m~ L ~ 



m 2 r ss 6M? w 12ml, . V CCB < V SM in this case if 

R 2 



1 M 2 h 2 v 2 ^ ( 1 - x 



Ak 2 m'jj^ 2 



< 



2x - 1 



Axil — x) 



(B17) 
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The r.h.s. has to be positive, i.e., x — > 0.5 + or mF t ole < 150 GeV. Slightly increasing rr? t ole 
would reverse the inequality. 

Lastly, there are two caveats. The first is that we used in this section the (one-loop 
improved) tree-level potential. However, following Ref. [49] we perform the calculations at 
the ^-scalar scale to eliminate large loop corrections. Secondly, a GCCB minimum may be 
"safe" if separated from the local standard-model minimum by a tunneling time greater than 
the age of the universe pi. Such considerations are beyond the scope of our present work 



and would also require a consideration of finite-temperature effects. 



APPENDIX C: CALCULATION OF THE LOOP-INDUCED MASS 



1. The EPM: Run and diagonalize 



A straightforward way to calculate A^o is (a) to take the first and second derivatives of 
AV(Q) with respect to the neutral CP-even and CP-odd components of Hi at the minimum; 
(b) to absorb the latter in m^o; and (c) to calculate the correction to (jLT]): the effective 
potential method (EPM). Like Ref. fPfl , we will follow Ref. |45|1 , however, we add to their 
expressions "D-term" [K| and Higgs-Higgsino and gauge-gaugino [^] contributions, which 
are typically —(1 — 2) and —(2 — 4) GeV, respectively. Note that (i) all the parameters 
are taken at the subtraction scale Q. The mass matrices are thus calculated at the scale 
Q and (ii) only then are diagonalized, i.e., a "run and diagonalize" algorithm. Our two 
choices in section [V| were Q = Mz and Q = 600 GeV ~ m^. (Large In ^gj- logarithms 
are multiplied by small couplings and do not invalidate the loop expansion, i.e., the EPM. 
This is equivalent to the statement above that Higgs-Higgsino contributions are small.) One 
corrects for threshold effects by the subtraction of one-loop leading logarithms (included in 
AV), and thus to that order the two choices are equivalent. 

The Q dependence of rriho(Q) comes about from (i) logarithmic wave-function de- 
pendence of i>i{Q) and (ii) two- loop implicit field- independent dependences of AV(Q) = 
AV[rrii(Q), QZi(Q), h t (Q)} (where rrii runs over the relevant mass parameters). Both are of- 
ten underestimated. The importance of (i) is obvious, e.g., from the m\jM\ factor. The 
field- independent dependence (ii) was studied in Ref. [pEfl , where it was shown to be of 
the order of two-loop (next-to-leading) logarithms. One can instead examine directly the 
expression for A^o, e.g., the leading logarithm in ([14]) evolves with scale (neglecting A# t and 
multiplicative factors) roughly as 



In 



m?(<2 2 ) 



o7T 



m 2 AQ 2 ) 



+ 



128vr 4 



Kh 2 t 



+ 2 



3 Us^? r in 



(CI) 



where k — 1, 2 for tL,R, respectively, and we took m t (Q) w rn t (Qo), etc., m~ th ~ m^, and 
neglected all terms aside from h t and g 2 s (= 47ra s ) ones. (Mg is the gluino mass.) This 
is of course a rough estimate, as, in practice all the parameters scale with Q. However, it 
illustrates qualitatively the presence of residual two-loop logarithms which are enhanced by 
the large couplings and masses when evaluating A^o. Note that typically the QCD term 
wins and the correction is positive near Q = Mz- 

Though formally of two-loop order, the residual Q dependence of the EPM m h o is (in 
our case) of order 10 — 20%. The ambiguity in m h o is reduced for m£ ^> 0. Also, for values 
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of tan/5 outside the allowed region Ag t is diminished (e.g., in the notation of Ref. jE| only 
A22 is important), h t < 1 (which also diminishes the wave-function renormalization) , and 
thus, there is a smaller ambiguity not only in myfi but also in A^o . Nevertheless, one should 
be aware of such ambiguities, which are generically present in a calculation of that sort. 



2. The RGM: Diagonalize and run 

One could alternatively correct for thresholds by explicitly decoupling a particle below 
its threshold and redefining the effective field theory: the renormalization group method 
(RGM). Following Ref. |]36| , |47| , we diagonalize the Higgs mass matrices at scale Q = m^o, 



integrate out the heavy Higgs doublet, and define a one Higgs doublet model (1HDM) 
effective theory with a quartic coupling |/i 4 . Note that m^o 3> mj, so only the heavy Higgs 
doublet is integrated out at Q = m^o, and other heavy particles will be decoupled at lower 
thresholds. The boundary condition for A(m^o) is determined by the MSSM Higgs couplings 
at that scale and requires knowledge of tan/^m^o) (i.e., wave-function renormalization). 
Mixing effects are accounted for by a complicated set of matching conditions [HS|,E7[. A is 



then evolved down to Mz, where m^o = y/Xv = 2\f\Mw I g 2 \ i.e., a "diagonalize and run" 
algorithm. 

One can solve the RGE's iteratively to get a semi-analytic approximation p6|j4T||. This is 



not an exact procedure, and at one- loop is accurate to leading logarithms (i.e., the accuracy 
of the EPM). Alternatively, if the RGM calculation is done by exact (or numerical) inte- 
gration, then the leading logarithms are summed to all orders. In agreement with Ref. |36[ 



we find that \{Mz) (and thus m^o) is diminished when integrating the RGE's numerically. 
However, the effect of the numerical integration (Fig. |5|) is more important than anticipated 
in Ref. p6[ since A is heavier than the mass-scale used there and because of the large cou- 



plings. The difference between the numerical treatment and the approximate formula is of 
order of two-loop (leading) logarithms. The summation of the leading logarithms modifies 
the scale dependence of the calculated mass. 

A different scale dependence is via the details of the decoupling of heavy states. To ease 
the calculation, we do not diagonalize scalar-quark mass-squared matrices at an intermediate 
scale. Instead, we put in the appropriate matching conditions f48| , |47|j by hand. These 



conditions are derived by expanding the EP assuming \m t fi\ ^ m^- , etc., which still holds 
in our case [m^m^o) < m t (Mz)]- Nevertheless, we repeated the 1HDM RGM calculation 
but using the EPM with Q « m^o to determine 77^0 (ra^o), and thus [using ^(m^o)] the 
boundary condition for A(m^o). The agreement with the calculation using the boundary 
conditions of Ref. pE5|,pET| is good. 



3. Validity of the comparison between the EPM and RGM 

Whether "running and diagonalizing" or vice versa leads to a slightly different answer. 



However, it was already noted in Ref. |36| that the difference between the two algorithms 
is small and of higher order, and can be ignored for the purpose of our discussion here. 
However, the inclusion of the appropriate RGM matching conditions is essential for the 
comparison between the two methods. 
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More importantly, the RGM approximation of 1HDM breaks down if is not sufficiently 
large. (In such cases the more conventional 2HDM RGM is appropriate.) Too small scalar- 
quark (diagonal) masses break down the expansion of Ref. 0,0] and too large scalar 
threshold-corrections can break down the EPM approximation. Indeed, one has to ensure 
that the scalar quarks are not heavier than about 1 TeV in order not to invalidate the loop 
expansion which underlies the EP, as these large logarithms would be multiplied by large 
couplings^ (see, for example, Ref. |57| ). We checked that the points in the parameter space 
where the RGM and EPM are in sharp disagreement (e.g., in Figs. [5] and ^) correspond to 
either one of these cases. 



4. Two-loop calculations 

In this work we do not carry out a two- loop calculation of A^o. Such a calculation is 
an elaborate task and attention should be paid to various issues, e.g., wave-function renor- 
malization and the coupled running of the masses and couplings; pole masses vs. running 
masses; the multi-scale structure of the low-energy theory; the correct choice of matching 
conditions (logarithmic and non-logarithmic) given a choice of decoupling scales; one-loop 
finite pieces (not included in the EP) that may become important in certain regions of the 
parameter space; in addition to correct counting of powers of the couplings. In fact, an 
hybrid algorithm between the EPM and RGM is required. 

We received recently two interesting papers where (partial) two-loop analyses were pre- 



sented |p8| , |59| . In this section we compare those results with the residual two- loop terms we 
identified above. We will take a somewhat critical point of view emphasizing that though the 
above studies (as well as the discussion here) pave the way to understanding A^o beyond 
one-loop order, they are not complete. We conclude that the one-loop leading-logarithm 
result is most probably an upper bound on the two- loop result, but we doubt whether any 
stronger conclusions can be drawn at the present. 

Above, we presented the leading- logarithm result. In the EPM this is given by / \v one ^ loop 
and in the RGM one solves the RGE's to that order. We found no significant difference 
between the two methods. Next, we proceeded to calculate beyond the leading order and 
integrated the RGM RGE's numerically, summing the leading logarithms to all orders in 
perturbation theory, and, in particular, to two-loop order. We then found a ~ —10 — 20% 
correction to A^o. On the other hand, we identified in the EPM expressions next-to- leading 
two-loop logarithmic terms. Those terms come about from the scale dependence of the 
running parameters, and introduce a similar (but positive) correction to A h o. Also, we 
noted the important (~ 10%) effect of ml° le / ra t unnin9 ~ [1 + (4a s /37r)] [e.g., ( m ™ nnm9 ) 4 pa 
0.8(mj° Ze ) 4 ]. It is then suggestive that the two-loop leading logarithms and m t corrections 
would diminish the one- loop (leading-logarithm) result, while the next-to- leading logarithms 
would partially counter-balance that effect. Thus, we reach the conclusion that the one- loop 
leading-logarithm prediction for m h o (which contain some next-to-leading logarithms) is an 
upper bound on the two-loop prediction. Any stronger conclusion would require a detailed 
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This requirement is consistent with our choice of ranges for the soft parameters. 
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consideration of the issues listed above. Furthermore, the significance of a detailed two- 
loop study in the context of SGUT's is not clear unless complemented by similar studies of 
threshold effects near the high-scale boundary. 

Let us proceed with the issue of the orders in perturbation theory. It was already pointed 
out by Chankowski J5?| that the inclusion of ^y one ~ loop in (0) with only one-loop RGE's for 
the (RG-improved) potential parameters is an inconsistent procedure when counting powers 
of couplings. These issues have recently received great attention in the context of theories 
with a single scalar field (i.e., 4 or the SM) |60,6l|. It was explicitly shown that a RG- 



improved L-loop effective potential is consistent when the parameters are calculated using 
(L+l)-loop RGE's 0. 

The observations of Ref. H60| , |6l[| were recently applied to the MSSM by Kodaira et al. 



who extended a previous work of Espinosa and Quiros ||62|| . Indeed, they find that the 
two- loop RGM with one-loop EP slightly increase the one- loop RGM result. However, they 
impose the naive assumption of only one relevant scale {Msusy) below which the SM RGE's 
are in effect (and the EP contain only a m t correction). Thus, the MSSM parameters are 
not treated on the same footing as the SM ones, and one expects some modification of their 
numerical results. However, we note that the correction is positive: recall that the difference 
between the two RGM approximations is the next-to-leading logarithms, which we found to 
be positive by examining the scale dependence of the result. 

Hempfling and Hoang [59] explicitly calculate the two-loop RGM correction (using the 
methods of Ref. f47fl), finding 
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9 2 s ~ 5ht 



(C2) 



which roughly agrees with our Eq. (|CT|). 

Those authors go further to derive a diagramatic two-loop result (in a simplified model) 
which diminishes A^o below its two-loop RGM prediction. They then show that, e.g., if 
the difference between the pole and running t-quark mass is accounted for in the RGM 
calculation by appropriate matching conditions, the two methods agree. This again stresses 
the need for adequate choice of matching conditions. However, it is not clear that there are 
no other effects of that order of magnitude that can slightly increase A h o. [Also, we showed 
that scalar mixings which are neglected in the above two-loop calculations and in our Eq. 
(PT|), are n °t apriori negligible.] 

In short, all analyses support the conclusion that two- loop logarithms will diminish the 
one-loop leading logarithm result (in particular, if one properly distinguishes the t-quark 
pole mass from the running mass). However, it is not clear that the prediction of a complete 
two- loop analysis lies as low as suggested in Ref. |59| (or by our Fig. |10|). (If indeed this is 
the case our upper bounds will be significantly strengthened.) 



APPENDIX D: EXTENDED HIGGS SECTORS 



If there is an extended heavy Higgs sector (as is assumed by some authors when trying to 
explain the light fermion spectrum) with large representations coupling to the third family 
(i.e., there is no flavor symmetry that forbids such couplings) then (fj) does not hold and our 



24 



constraints are evaded. Alternatively, some models extend the light Higgs sector by adding a 
SM singlet superfield, S, to the spectrum. This is the only addition to the Higgs sector which 
is consistent (i.e., without fine-tuned cancellations) with coupling constant unification. 

The new singlet (s) would mix with the other neutral Higgs bosons and would thus con- 
tribute to m^ ; e.g., fl63[j , 

2 

mlo 2 < M|cos 2 2/? + A 2 y sin 2 2/?, (Dl) 

where A s (« a ) is a new Yukawa coupling in the superpotential: X S H 1 H 2 S (k s S 3 ). (We do 
not distinguish here our notation of superfields from that of their scalar components.) The 
MSSM limit with m^ = is recovered if A s — > 0, which is indeed the case if h t ~ h{ lxed ||64|| . 
However, the new Yukawa couplings (i.e., X s ^ 0) enable one to realize (0) for somewhat 
smaller values of ht < h{ lxed (and tan/3 > 1). The shift in the allowed region is shown in 



Fig. |TT] for A s = 0.5 and k s = 0. Thus, m^ can be large and m^o can be heavy ~ 150 GeV 



T9| , |20|j6^| when loop corrections are added. (See also Ref. 
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FIGURES 



FIG. 1. The i-quark pole mass (in GeV) - tan/3 plane is divided into five different regions. 
Two areas (low- and high-tan (3 branches) are consistent with perturbative two- Yukawa (b and r) 
unification and with an upper bound of ^|| GeV for the 6-quark current mass. The 0.85 factor takes 
into account a ~ ±15% theoretical uncertainty in the numerical calculation due to the matching 
conditions. Between the two branches the 6-quark mass is too high. For a too low (high) tan /?, 
the t-quark (6-quark) Yukawa coupling diverges. The region where all three (third- family) Yukawa 
couplings unify (dashed line) intersects the allowed high-tan (3 branch. [The three- Yukawa region 
is calculated without imposing any constraints on mj, and is assigned a ~ 5% uncertainty (shaded 
area) from corrections to the h t /hb ratio.] The i-quark mass range suggested by the electroweak 
data and the tan (3 = 1 line are indicated (dotted lines) for comparison. 

FIG. 2. The prediction for the Higgsino mass parameter /j, (after absorbing AV in the redefined 
tree-level parameters) as a function of (a) the t-quark pole mass and of (b) tan/3. (The t-quark 
mass and tan (3 are strongly correlated.) Filled squares indicate (probably unacceptable) points for 
which the SM minimum is only a local minimum. All masses are in GeV 

FIG. 3. A scatter plot of the Higgs boson mass (calculated in the EPM with the subtraction 
scale Q at the Z-pole) vs. the t-quark pole mass. The o indicate points which have no (potentially 
dangerous) negative-energy color and/or charge breaking minimum. The o have such a minimum 
but it is a (safe) local one, i.e., it lies above the standard- model minimum. Points indicated by a 
filled diamond have an unacceptable CCB global minimum. The Z mass is indicated for comparison 
(dashed line). All masses are in GeV. 

FIG. 4. Same as in Fig. 3 except the mass is calculated using the RGM (numerically). 

FIG. 5. The EPM (with the subtraction scale Q at the Z-pole) calculation of the Higgs 
boson mass is compared with the mass calculated in the RGM but using the Haber and Hempfling 
leading logarithm formula. All masses are in GeV. Points corresponding to global GCCB minima 
(indicated by a filled diamond in Figs. 3 and 4) are omitted. 

FIG. 6. Same as in Fig. 5, except Q = 600 GeV and the RGM calculation is carried out 
numerically. 

FIG. 7. The Higgs boson mass (in GeV) distribution in a sample monte-carlo calculation 
(using the EPM with the subtraction scale Q at the Z-pole) for the t-quark pole mass in the range 
(a) [155, 165], (b) [165, 175], (c) [175, 185] GeV. 
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FIG. 8. Same as in Fig. 7 except for the omission of points which correspond to a global color 
and/or charge breaking minimum (marked by a filled diamond in Fig. 3). For comparison, the 
shaded areas indicate points which are consistent with the constraint (13) (marked by a circle in 
Fig. 3). The difference between the two distributions is points with only a local CCB minimum 
unnecessarily excluded by (13) (marked by open diamonds in Fig. 3). 

FIG. 9. Same as Fig. 7 except for substituting the t-quark running (rather than pole) mass in 
the mass expressions for the scalars. 

FIG. 10. Same as in Fig. 7a-b except using the RGM (with numerical integration). The 
distributions omitting points which correspond to a global color and/or charge breaking minimum 
(marked by a filled diamond in Fig. 4) (lightly shaded areas) are also indicated, and should be 
compared with the total (shaded and unshaded) areas in Fig. 8. (For a t-quark heavier than 
about 180 GeV the 1HDM assumption that underlies the calculation is not always accurate - see 
Appendix C.) 

FIG. 11. The low tan/3 branch in the presence of an additional Higgs field with X s = 0.5 and 
k s = 0. The allowed area when there is no singlet field (A s = 0) is shown for comparison (shaded 
region) . the t-quark pole mass (given in GeV) range suggested by precision data and tan [3 = 1 are 
indicated (dotted lines). 
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TABLES 



TABLE I. The allowed regions in the tan/3 



m 



pole 



plane. Relations (1) and (2) are also 



consistent with the parameter space described by case (2) except that ht > h^. rriq stands for 
a typical scalar-quark mass. Finite corrections oc tan/3Afi /u/mg could modify the allowed m^ ole 
range in case (2). 



case (1) 



case (2) 



GUT relations 
The allowed m\ ole range (GeV) 
The allowed tan (3 range 
Yukawa couplings 
The Higgs potential (approximate) symmetry 
The vacuum (approximate) symmetry 
SU(2) x U(l) breaking "messenger" 
The large mass parameter 



r 
h° 



A h o enhancement 
A , H+ 



© and (|) 

140 190 

1 - 2 

ft* » 1 » fc 6 
SU{2) L x SU(2) R 
SU(2) L+R 
ml 

Higgsino mass parameter 

left-right mixings 

mjjo ~ m^o 
heavy 
m A o « V2\n\ » 2 TeV 



© and (|) 

~ 180 ± 15 

50 ± 10 
h t ~ h ~ 1 

4 x 4 

scalar-quark masses 



large m< and the heavy scalars 
massive pseudo-Goldstone bosons 
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Abstract 

The SU(5) unification-scale relation hf, = h T between the 6-quark and r- 
lepton Yukawa couplings severely constrains tan (3 and the i-quark mass (even 
more so if h t = hf, = h T holds) in supersymmetric models. We examine the 
implications of these constraints for the Higgs sector assuming universal soft 
breaking terms, and emphasize that both of these relations impose unique 
characteristics in terms of symmetries and of the spectrum. We further study 
the tan/3 ~ 1 scenario, which is suggested by hf, = h T , and, in particular, the 
loop-induced mass of the light Higgs boson. We compare the effective poten- 
tial and renormalization group methods and stress the two-loop ambiguities 
in the calculation of the mass. These and a large enhancement to the loop cor- 
rection due to i-scalar left-right mixing considerably weaken the upper bound 
on the mass of the light Higgs boson that has been reported. Nevertheless, we 
find that for this scenario the Higgs boson is probably lighter than 110 GeV, 
and typically lighter than 100 GeV. Thus, it is in the mass range that may 
be relevant for LEPII. Our numerical results are presented in a self-contained 
manner in section V. In separate appendices we discuss the global symmetries 
of the Higgs potential, the issue of false (color-breaking) vacua, which may be 
important for tan (3 ~ 1, two-loop calculations, and the effect of an additional 
Higgs singlet. We show that the approximate constraints that are often used 
to eliminate color-breaking vacua are not always relevant. 

PACS numbers: 12. 10. Dm, 11.30.Pb, 14.80.Gt, 14.80.Ly 
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I. INTRODUCTION 



Supersymmetric grand-unified theories [1,2] (SGUT), though far from being proved, are 
an attractive framework on both observational and theoretical grounds. (Their relation 
to string theory remains at present poorly understood.) In the minimal models, which 
we assume hereafter, there is a grand desert between the weak and the unification scales 
(M z ^ Q ^ TeV and M G 10 16 - 10 17 GeV, respectively). At the unification point M G 

ai {M G ) = a 2 (M G ) = a 3 (M G ) = a G (1) 

(where a.\ = |cty and a G 0.04), and, depending on the grand unifying group representa- 
tions of the Higgs superfields that couple to matter, certain relations between the standard 
model Yukawa couplings hold. In the minimal models based on SU(5), SO (10), and E 6 
these representations are the fundamental ones, and one obtains [3] 

h T (M G ) = h b (M G ). (2) 

For some models based on SO (10) (or larger groups) one has the stronger prediction 

h T (M G ) = h b (M G ) = h t (M G ), (3) 

but this only holds if one makes the additional assumption that the masses are generated by a 
single complex Higgs 10-plet. (h Tj b,t are the SM Yukawa couplings of the r, 6, t-superfields, 
respectively.) One usually assumes that some perturbation modifies the coupling or the 
masses of the two light families where, in principle, similar relations should, but do not, 
hold 1 . (See also Appendix D.) Also, such relations often do not hold in superstring unified 
models. 

At the weak-scale the particle content is that of the minimal supersymmetric standard 
model (MSSM), i.e., that of a two Higgs doublet model but with each bosonic (fermionic) 
degree of freedom complemented by a fermionic (bosonic) one. In particular, the Higgs 
sector contains three Goldstone bosons and five physical degrees of freedom - two CP even 
(h° and H°) and one CP odd (A ) neutral, and one complex charged (H + ) Higgs bosons. 
The symmetries and spectrum of the Higgs scalars are strongly affected by the above GUT 
assumptions, and are the subject of our present study. However, a review of previous results 
in SGUT's, which establish the base for our present work, is appropriate. 

The consistency of (1) with the data was pointed out recently in Ref. [4] and that of 
(2) and (3) in Ref. [5]. Relation (1), however, holds only if there is exactly one pair of 
Higgs doublets (but any number of SM singlets is allowed). We previously investigated the 
status of coupling constant unification in great detail [6] where we considered non-trivial 
matching conditions 2 that perturb relation (1) in a model dependent way. We concluded 



1 Since the Yukawa couplings corresponding to the third family are much larger than those of the 
two light families, relation (2) [or (3)] will be only slightly affected. 

2 Heavy degrees of freedom (which are model dependent) are integrated out of the effective theory 
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that no significant constraints on the MSSM parameters can be deduced just by that relation. 
However, we find 3 that (1) predicts for the strong coupling [using a(Mz)~ 1 = 127.9 ±0.1 [7] 
and s 2 (M z ) = 0.2324 ± 0.0003 (for m p t ole = 143 GeV) for the (MS) fine-structure constant 
and weak angle, respectively] 

a s (M z ) = 0.125 ± 0.001 ± 0.008 + H a , + 3.2 x 10- 7 GeV- 2 [(mr' e ) 2 - (143 GeV) 2 ] , (4) 

where the ±0.008 (±0.001) uncertainty in (4) is theoretical (due to the input parameter error 
bars). The former is due to the unknown values of the matching conditions (i.e., threshold 
and nonrenormalizable-operator effects) mentioned above. The function H as is a correction 
from Yukawa interactions. It is negative but negligible, unless some Yukawa couplings are 
<; 1, for which we obtain H as ~ —(0.001 — 0.003). The quadratic m p t ole dependence is 
induced by the correlation between the weak angle extracted from the data and the t-quark 
pole mass, m pole . Thus, typically one predicts a s (Mz) ^ 0.12 with a significant dependence 
on m pole . This observation is of importance when simultaneously considering (1) and (2) [or 
(3)] to predict the 6-quark (MS running) mass mj [and/or the t-quark mass, if using (3)] in 
terms of the r-lepton mass. We find [using, e.g., Eq. (4) but with a more careful treatment 
of theoretical uncertainties to the 6-quark mass prediction] that only a small region in the 
tan (3 — m pole plane is allowed by the experimental range of mj [8], where tan (3 = v U pl^down 
is the ratio of the two Higgs doublet expectation values. This holds even when non-trivial 
matching conditions at both scales are considered. 

It is customary to assume that the MSSM spectrum is given to a good approximation 
near Mq by a small number of universal supersymmetry-breaking soft parameters; i.e., a 
common scalar mass m , trilinear and bilinear (dimension-one) couplings A and B 0} and a 
common gaugino mass Mi; in addition to a supersymmetric Higgs mass (i.e., the Higgsino 
mass), fi [2]. We will adopt these assumptions below. (For an alternative scenario, see, 
for example, Ref. [9].) However, our conclusions so far are largely (i.e., except for the the 
details of the matching conditions) independent of those assumptions and hold in any SGUT 
scenario with an approximate grand-desert in which the MSSM /3-functions are valid. The 
universality assumption just strengthens the constraints on the tan/? — m pole plane; i.e., 
assuming universal initial conditions and requiring that the weak-scale MSSM (one-loop 
improved) scalar potential is consistent with a broken SU{2) X U(l) symmetry, we have to 
further constrain tan/? > 1 and h t > h b . We then find two allowed regions, 

1. tan (3 = 1 + A/?, 140 < m pole (< 200) GeV, h t « 1 > h b ; 

2. tan f3 « 50 ± 10, 170 ± 10 < m pole (< 200) GeV, h t > h b « 1. 



at Mq and some modification to the naive boundary conditions in (1) and (2) [or (3)] is required, in 
general. The appropriate matching functions may also include corrections from nonrenormalizable 
operators since Mq is not too far from the Planck scale. Similarly, integrating out the MSSM new 
particle spectrum will correct the boundary conditions at Mz- 

3 Eq. (4) updates the results in [6] for the effect of more recent precision electroweak data. 
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The first of these solutions is consistent only with hf, = h T , and has A^ < f for m P ° e < 190 
GeV. A part of the region which is described by the second solution is consistent with 
the stronger prediction h T = hf, = h t . A third solution with 3 < tan/? < 40 exists for 
m p t ole 215 ±10 GeV, which is the upper bound on m pole if the model is to stay perturbative 
up to Mq. However, such large values for m t are inconsistent with precision electroweak 
data. The allowed regions are illustrated in Fig. 1. 

These results are in good agreement with those of Ref. [10-14], with any small differences 
due to different treatments of matching conditions and of the coupling constant unification 
condition. Arbitrarily relaxing both condition (1) and (2) independently [e.g., a s (Mz) = 
0.11 and hb(Ma) = 0.85/j t (Mg)] could invalidate our conclusions. However, the corrections 
to (1) and (2) are strongly correlated, and when treated as such the strong constraints 
on the tan/? — m pole plane hold 4 . In particular, the m;,/m T ratio is relatively insensitive 
to the details of the decoupling of the low-scale spectrum [8], and our conclusions would 
not be significantly modified if those thresholds were treated differently than in our (semi- 
analytic) approximation (where we used six effective mass parameters). Let us stress that 
any treatment of theoretical uncertainties would be lacking unless the potentially important 
corrections at the high scale are accounted for in one form or another. These issues are 
summarized in Ref. [15,16], where the second reference updates some of the results of Ref. 
[6,8]. Here, relevant results are again updated using a preliminary analysis of most recent 
precision data. 

A two-parameter fit [allowing the light (SM-like) Higgs boson mass, m h o } to vary from 
50 — 150 GeV with a central value ra^o = Mz] to all Z, W, and neutral-current data yields 

s 2 (M z ) = 0.2324 ± 0.0003 - 0.92 x 10" 7 GeV" 2 [(mr' e ) 2 - (143 GeV) 2 ], (5) 

mr' e = 143i 19 7 + 12.51n^GeV, (6) 

which in turn implies tan/? ^ 1.3 (at one s.d.). Eqs. (5) and (6) update Ref. [17] where 
m pole = 134^28 + 12.5 In GeV was given (which also implied tan /? ^ 1.3). 

It was pointed out by Barger et al. [12] that if indeed m pole < 160 GeV (as suggested by 
the range given in Ref. [17] and here) then ra^o [in case (1)] is determined by the magnitude 
of the loop corrections, A^o, and is < 85 GeV. If so, the 6-quark mass prediction in minimal 
SU(5) (or similar) SGUT constrains the SM Higgs boson mass more significantly than the 
general MSSM (triviality) upper bound of ~ 130 GeV [18-20]. We agree qualitatively with 
this conclusion, but argue that two-loop ambiguities and important left-right t-scalar mixing 
effects, in addition to the higher m pole range, weaken the bound to m h o < 110 GeV. 

Below, we will establish the basis for and then pursue further numerical studies of the 
prediction for ra^o and its upper bound in the tan /? 1 scenario. Rather than approximating 
the one-loop correction we will calculate ra^o for a given point in the parameter space using 
the supersymmetric spectrum (calculated numerically), and use monte-carlo routines to 
study the upper bound. In section II we briefly review the MSSM Higgs potential and its 



4 That is unless one assumes a conspiracy of various corrections that are otherwise independent. 
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minimization conditions. We specialize the discussion to the minimal SGUT scenarios, cases 
(1) and (2), and emphasize that the Higgs potential exhibits approximate global symmetries 
in these cases, on which we elaborate in Appendix A. We concentrate thereafter on the 
tan f3 ~ 1 scenario. Some general features are described in section III. In particular, the 
large Higgsino mass parameter can generate false vacuum solutions with broken color and 
charge. We discuss that issue in greater detail in Appendix B. The calculation of the 
light Higgs boson mass in both the effective potential and renormalization group methods 
is discussed in section IV, where we use the comparison between the two methods to study 
the two-loop ambiguities in the calculation. The details of the two methods, as well as 
two-loop calculations, are further discussed in Appendix C. Section V is reserved for further 
numerical studies of the Higgs boson mass, and we show that because of two-loop ambiguities 
and II — Ir mixing m^o can easily exceed 85 GeV but, unless m p t ole is much heavier than 
the range suggested by (6), it remains in the range that may be relevant for LEPII. There 
is even some possibility that it is in the range still relevant for LEPI. We also discuss the 
relation between t-scalar mixing and color breaking. In particular, parameter ranges which 
correspond to large mixing (and large m h o ) often have an unacceptable color-breaking global 
minimum, but not always. The discussion in section V is to a great extent independent of 
the other sections and can be read on its own. We will summarize our conclusions in section 
VI, where we also examine the implications for m p t ole . In Appendix D we point out simple 
extensions of either the MSSM or the GUT in which our (constrained) analysis does not 
apply. 

Below, we point out that (negative-energy) color-breaking local minima of the full scalar 
potential are generic. Approximate analytic constraints that are often used in the literature 
attempt to eliminate not only global but also local color-breaking minima, and are thus too 
strong. Furthermore, they are designed to eliminate only certain types of color-breaking 
minima, and are thus too weak. This is explicitly demonstrated in Appendix B, where we 
also compare analytic and numerical treatments of the problem. 

We follow Ref. [8] in calculating the couplings (and the unification point), and Ref. 
[21] in treating the one-loop effective potential correction AV, including contributions from 
all sectors. The boundary conditions at Mq for < m < 500 GeV, |A | < 3m , and 
50 < Mi < 350 GeV, are picked at random. m pole values are picked at random but with a 

gaussian distribution defined by (6). For a given mf°' e , tan ft(Mz) is picked at random within 
its allowed region [case (I)]. Values of mf e < 155 GeV are discarded so that tan f3 ^ I.I 
and the divergent limit is not reached (see below). We then solve iteratively for //q (and 
B ). [The relevant renormalization group equations (RGE's) have been given by various 
authors [22].] The sign of /i , which is a RG invariant, is picked at random as well. We 
explicitly verify that the potential is bounded from below; SU{2) X U(l) is properly broken; 
the lightest supersymmetric particle (LSP) is neutral (in practice, it is a neutralino); and 
that all the squared masses which correspond to physical scalars are positive before a point 
is accepted. We then find, in general, that we are in agreement with lower bounds coming 
from direct searches, (e.g., values of Mi < 100 GeV, which would imply a too light chargino, 
are already excluded in our case by a tachionic t-scalar.) We do not impose a lower bound 
on m h o } which is the subject of our investigation. Our treatment of color and/or charge 
breaking (CCB) minima is discussed in Appendix B. 
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II. THE WEAK-SCALE HIGGS SECTOR 



The Higgs part of the MSSM (weak-scale) scalar potential reads [23,18] 



V{H U H 2 ) = (m 2 Hi + + (m 2 H2 + ^ 2 )\H : 



2 



\ MSSM 

+£//(#i# 2 + h.c.) + (\H 2 \ 2 - m 2 ) 2 + AV, (7) 

where m 2 H , rn 2 H 2 ^ an d B (//) are the soft (supersymmetric) mass parameters renormalized 

down to the weak scale, ra\ = B\i < 0, A MSSM = 92+ £ 91 ; and we suppress SU(2) indices. 
The one- loop correction [24,18] AV = AV one -' oop (which, in fact, is a threshold correction 
to the one-loop improved tree-level potential) can be absorbed to a good approximation in 
redefinitions of the tree-level parameters [21,25]. 

A broken SU{2) X U(l) (along with the constraint m 2 H -\-rn 2 H2 +2/i 2 > 2|m 2 | from vacuum 
stability) requires [23,18] 

K fi +// 2 )K f2 +// 2 )<|m!| 2 , (8) 
and the minimization conditions then give [23,18] 

9 m 2 rr — m 2 rr tan 2 8 1 9 . . 

2 = _Jh H2 1± _ _ M 2 9a 

P tan 2 /3-l 2 z ' y J 



m 2 . 



-isin2/3[m^ i +m^ 2 +2 /U 2 ]. (9b) 



For tan/? — > 1 one has h t ^> hf, } and hence cannot have m 2 H = rn 2 H2 (assuming universal 
initial conditions). Thus, \fi\ — > oo in that limit [case (1)], and the SU{2) X U(l) breaking is 
driven by the Bfi term. In practice, we will stay away from the divergent limit in case (1) by 
taking tan 8 ^ 1.1 [21]. We do not expect significantly different results for 1 < tan f3 < 1.1; 
i.e., the divergence is either stabilized by model dependent finite-loop corrections or would 
exclude that region (see also below). For tan f3 — > oo [case (2)] one has Bfi — > so that the 
symmetry breaking is driven by m 2 H < 0. 

The tan f3 — > 1 case corresponds to an approximate SU (2)l + r custodial symmetry of 
the vacuum [26,27], which we further discuss in Appendix A. The symmetry is broken at 
the loop level so that one expects tan f3 slightly above unity, in agreement with our cut 
tan f3 <; 1.1. As a result of the symmetry, the CP-even Higgs mass matrix becomes 

M 2 « tJ 2 x ( _ X X ^ , (10) 

and it has a massless tree-level eigenvalue, mf ~ 0. This is, of course, a well known result 
of the tree-level formula [23,18] 



T 2 1 



m 2 A0 +M 2 Z - ^J(m 2 A0 + M 2 ) 2 - im 2 A0 M 2 cos 2 2/3 



'ID 
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when taking /? = ^. The mass is then determined by the loop correction ra^o sa A^o oc h t m t 
and is further studied in sections IV and V. On the other hand, in case (2) /? — > -| and we have 
mjp ~ Mz (assuming m^o > Mz)- When adding the loop correction 5 ra^o < y/2Mz ~ 130 
GeV. The heavier CP-even Higgs boson mass eigenvalue equals approximately \/2|/w| in case 
(1) and m H o sa m A o in case (2). The loop corrections are less relevant here as typically 
m 2 H o ^> A 2 H o (in particular, when \fi\ — > oo). 

The custodial symmetry (or the large // parameter) dictates in case (1) a degeneracy 

y/2\ii\. [The t ree-level corrections to that relation are of order 
(Mw, z/'m-A ) 2 -] That is, at a scale Q sa \/2|/w| ~ 2 TeV the heavy Higgs doublet H is 
decoupled, and the effective field theory below that scale has only one massless (mf ~ 0) 
Coleman- Weinberg SM-like {y^o = v) Higgs doublet, h. Unlike in our case, in general 
Mz/m^o is not guaranteed to be small, and the decoupling scale Q is much lower, typically 
below the scalar-quark thresholds. 

The Higgs sector in case (2) exhibits an approximate O4 X O4 symmetry (see Appendix 
A). This is a special case of the m\ — > case discussed in Ref. [28]. To obtain m\ sa 
(starting with universal initial conditions and hj, sa h t ) one would require large Mi and 
|/io|, which enable one to obtain the desired limit by adjusting large cancellations [29-31]. 
Once expectation values are acquired, A and H + are massive pseudo-Goldstone bosons 
of the broken O4 symmetry. However, they are typically heavy (for v& own = the O4 
symmetry is restored) but lighter than the scalar quarks in this scenario (whose 
^jml + 6M? w 1 TeV). 

Large tan/? [case (2)] solutions are discussed in Ref. [29,30,13,32,31]. They predict 
m p t ole sa 180 ± 15 GeV [8,13], heavy scalars, and contain potentially troublesome loop cor- 
rections, e.g., to m&, b — > 37 [13,32]; and, in the absence of light scalars to counterbalance 
large m t contributions, also to Z — > bb (see Ref. [33]). The former can eliminate the motiva- 
tion for any Yukawa coupling analysis. A different approach was recently presented in Ref. 
[31] where those large finite-loop corrections were used to construct models consistent with 
universality, a broken SU{2) X £7(1), and h t = hj, = h T . Finite superpartner (~ 1 TeV) loops 
are adjusted to diminish mj by the right amount (~ 1 GeV), and mj is strongly dependent 
on the values of the soft parameters, a possible but undesirable situation. Much lower values 
of m p t ole are then preferred. Diminishing m t may rectify the situation with Z — > bb, but the 
prediction for the b — > 37 branching ratio can be significantly modified by the same super- 
partner loops [13]. In contrast, the tan/? ~ I [case (1)] solution constrains only m pole > 140 
GeV (mr e > 155 GeV when imposing the tan /? > 1.1 cut), is in better agreement with (6), 
and predicts a light Higgs boson (ra^o ^ 110 GeV). Cases (1) and (2) are summarized and 
compared in Table I. We hereafter study the tan /? sa I [case (1)] solution only. We assume 
that finite-loop corrections to mj are at the most a few percent which is, in general, the case 
away from the large tan /? limit, and thus do not alter our results. [The parameter space is 
sensitive to those corrections only in case (2).] 



5 A/ l o is enhanced in this case by the large nit an d the heavy i-scalars. 
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TABLE I. The allowed regions in the tan (3 



m 



pole 



plane. Relations (1) and (2) are also 



•niq stands for 

pole 



consistent with the parameter space described by case (2) except that h t > h 
a typical scalar-quark mass. Finite corrections oc tan (3 M i fi I ml could modify the allowed m p t 

2 

range in case (2). 



case (1) 



case (21 



GUT relations 
The allowed m pole range (GeV) 
The allowed tan (3 range 
Yukawa couplings 
The Higgs potential (approximate) symmetry 
The vacuum (approximate) symmetry 
SU(2) X U(l) breaking "messenger" 
The large mass parameter 



m h0 

A h o enhancement 
A , H+ 



(1) and (2) 

140 190 

1 - 2 
h t « 1 > h b 
SU(2) L X SU(2) R 
SU(2) L+R 
m\ 

Higgsino mass parameter 
m\ « 
left-right mixings 
m H o ~ m A o 
heavy 
m A o « y/2\fi\ « 2 TeV 



(1) and (3) 
~ 180 ± 15 

50 ± 10 
h t « h b « 1 
4 x 4 
3 x 3 

m h 

scalar-quark masses 
m T h0 « M z 
large m t and the heavy scalars 
m H o ~ m A o 
massive pseudo-Goldstone bosons 
m A o < mg 



III. THE tan/3 ^ 1 SCENARIO 

Our primary motivation to study the tan f3 ~ 1 solution to the MSSM is its consistency 
with condition (2). The large values of a s (Mz) implied by (1) [e.g., Eq. (4)] generate (too) 
large positive loop corrections to hb/h T} which need to be counterbalanced by large negative 
corrections oc h^, h\. The latter are sufficient only if some of the the Yukawa couplings are 
large, i.e., near their unitarity upper bounds, which is the case in the allowed regions. [Also, 
H as in (4) becomes non-negligible in that limit.] The large values of the Yukawa couplings 
can be understood in terms of a quasi-fixed point in the flow of the respective RGE's [34]. 
In particular, in case (1) h t is near its quasi-fixed point. The two issues, i.e., the consistency 
of tan f3 ~ 1 with (2) and the /j r RGE flow structure, are strongly related. 

The smallness of (or alternatively, the large \fi\) can be understood in terms of an 
approximate SU (2)l x SU (2)r symmetry in the Higgs potential, which, if it exists at some 
scale, will be only slightly broken at Mz, as is shown in Appendix A. The tan f3 ~ 1 
scenario was previously studied in Ref. [35], where it was referred to as "highest classical 
degeneracy": in the tan (3 = 1 limit condition (8) becomes an equality. [In practice, we find 
that {m 2 H + /« 2 )(?7^ 2 + n 2 )jra\ > 0.99 for < 1 and after the proper redefinitions.] It 
was subsequently studied in Ref. [36] in the context of a global MSSM (i.e., no high-scale 
assumptions), and its consistency with limits on m h o was shown. The tan (3 ~ 1 solution 
suggested by the minimal SGUT relation h b = h T revived the interest in that corner of 
parameter space. More recently Ref. [37-39], motivated by the quasi-fixed point prediction 
to m t , extended the discussion from the Higgs sector to the full parameter space. 

In Figs. 2a and b we show the prediction for function of m p t e and of tan /3, 

respectively. Typically \fi\ ~ 1 TeV, depending on A^ and on the soft parameter scale. It is 
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the large \fi\ parameter that dictates the characteristics of the scenario (and not only for the 
Higgs sector). Before further elaborating on the calculation of ra^o, let us briefly summarize 
some of the features that appear in our numerical studies, and which are of relevance for the 
discussion below. Though we agree on the characteristics with other authors, our emphasis 
and interpretation are different. 

f. The Higgsino (with mass ~ //) is decoupled from the (much lighter) gauginos. The 
bino and the wino are the lightest neutralino and chargino, respectively. (The former 
is the LSP which has a bino fraction near unity.) The heavy Higgsino decouples from 
the Z — >■ bb vertex, and thus, smaller values of m p t ole are favored 6 . However, m p t ole 
and tan/? are correlated, and as m pole grows |//| is diminished (see Fig. 2). Thus, in 
principle, Higgsino - t-scalar loops can still counterbalance large m t contributions to 
Z — >■ bb (see, for example, Ref. [33]) if the t-scalar is light enough. [A light t-scalar is, 
however, less likely if \fi\ is small (see below) and from the /^-parameter.] 

2. The t-scalar mass-squared matrix is (we take all parameters to be real) 

m~ m t (A t + fi/ tan f3) \ 



m t (A t + fi/ tan f3) m 



2 



A2) 



and similarly for the superpartners of the other fermions. The diagonal elements 
correspond to the left and right-handed t-scalar squared masses, ^, which consist 

of soft, F, and D terms. A t is the trilinear soft parameter preceeding the hiF^t^R 
term in the scalar potential, with A t (Ma) = A . The prediction for the diagonal terms 
has no unique charachteristics in our case, and the left and right-handed masses are 
in the ~ fOO GeV — f TeV range. However, the mixing term is typically large (unless 
the fermion ITlclSS IS ~ as is the case for the light families). In the limit \fi\ — > oo 
there is a 6 t = ^ left-right mixing so that one of the eigenstates would have a negative 
mass squared, and the limit cannot be fully realized. In practice, the soft terms and 
m t |/i| are often of the same order of magnitude, with a nearly degenerate mass-squared 
matrix. Thus, the requirement that (9) defines a physical minimum with no negative 
squared masses, corresponding to physical scalars, constrains //. Also, our previous 
comment regarding a light t-scalar and a small //-parameter is clear by observation of 
(12). It is interesting to note that by fine adjustments of the soft parameters [i.e., of 
the degeneracy in (12)] one may realize a scenario with a light (< 45 GeV) t-scalar 
which is nearly decoupled from the Z (e.g., sin# t 0.8). (See, for example, Ref. [41].) 
However, such scenarios are often associated with false vacua (see below and Appendix 
B). 

3. The absence of a negative t-scalar squared mass determines a (m and dependent) 
lower bound on Mi. In particular, "no-scale" (i.e., m = A = 0) SGUT models are 



3 The prediction for the Z —> bb branching ratio decreases with m pole , and its measured value is 
' 1.5 s.d. higher than the SM prediction with m P ° e ~ 140 GeV. The discrepancy grows with m P ° e 



[40]. 
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strongly constrained 8 . Also, requiring that the light eigenstate of the r-scalar mass- 
squared matrix (which carries charge) is not the LSP determines an upper bound on 
Mi in the no-scale limit (see also Kane et al. [20]). These point at a correlation between 
the bounds on the soft parameters. A different class of correlations which take effect 
in the h t — > h{ txed limit is discussed by Carena et al. [37]. Strong correlations (e.g., 
the m p t ole — tan (3 correlation 9 , the h t — > h{ txed ones, and those due to scalar mixings 
discussed here) do not allow elimination of the large- \fi\ solutions (that correspond to 
large left-right mixing enhancement to m^o - see below) on the basis of fine-tuning 
arguments. The naive notion of fine-tuning becomes obscure and misleading in the 
presence of strong correlations. (See also Ref. [37].) Rather, strong correlations make 
the scenario quite predictive. 

4. We can rewrite the mixing term in (12) as hfH^T where = v U p- T is strongly 
constrained by requiring that the full scalar potential does not have a color and/or 
charge breaking global (GCCB) minimum, which is a different issue than whether the 
extremum defined by (9) is indeed a minimum. The issue of false vacua is of particular 
interest in the tan /? — > 1 scenario as V(Hi } H 2 ) — > 0~ in that limit. Below, we will 
consider that subject in some detail. In particular, see Appendix B. We find that the 
absence of negative energy CCB minima implies the constraint, first described in Ref. 
[43], 

(|A t | + 5 H) 2 <2(m? i+ my, (13) 

where s = A t fi/\A t fi\, (\A t \ + s|/i|) 2 = T 2 (for tan/? = 1), and m 2 ^ ^ consist of only 
the soft terms. (Note the importance of the trilinear parameter and of the relative 
sign.) Our numerical studies imply that (13) is (to a good approximation) a sufficient 
condition (in that region of parameter space) but is not necessary (i.e., it is too re- 
strictive). For example, more than 80% of the points that are inconsistent with that 
constraint for m p t ole > 165 GeV correspond to CCB minima which are only local and 
are therefore safe. In Fig. 2 only filled squares correspond to (probably unacceptable) 
GCBB minima. Also observe in Fig. 2 that // < is preferred by CCB constraints. 
This is because typically A t ^> A for A < so that the contribution from a large 
and positive // cannot be canceled in (13) (for our range of A ). 

IV. THE LOOP-INDUCED MASS 

Let us now proceed to discuss the mass of the SM-like higgs boson. The upper bound 
on m^o is given, e.g., in Ref. [44], 



The strict no-scale assumption Bo = uiq = Aq = is not consistent here. See Appendix A. 
8 No-scale models with hf, = h T were recently discussed in Ref. [42] 

9 Note that even if we do not require hf, = h T , m pole and the lower bound on tan/? are correlated 
due to unitarity considerations. 
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where 



m\ 



< Ml 



cos' 



2/3 + 



3am^ 



4tts 2 1 



In 



mf 



+ A 6 



(14) 



raj 



+ 



9 

raf 



9 

m| 

12 



ra 

sin 2 2fl t 
4m? 



*2 



2m 



-mi 

*2 



In 



m? 

*2 



(15) 



and where mf are the eigenvalues of the t-scalar mass-squared matrix, 6 t is the mixing 
angle, and we have neglected other loop contributions. (They are included in our numerical 
calculations below.) The tree-level mass squared, ra^ 2 , and the loop correction, A^ , are 
bounded by the first and second terms on the r.h.s. of Eq. (14), respectively. In the absence 



of mixings Ag t = 0. For tan (3 — > 1 one obtains 



ra 



h° 



0, and thus m h o A 



h°- 



The upper bound on A^o scales as m\ and depends on the overall scale of the t-scalar 
mass, the amount of left-right mixing, and the separation |ra~ — mjj- The large mixing and 
separation for \fi\ — > oo enhance A^o significantly (typically by 10 — 20 GeV) compared to 
the first term alone (which was all that was included in Ref. [12]). Condition (2) (and h t 
unitarity) correlate m p t ole and tan/? in case (1) (see Ref. [8,16]), and A^ increases from ~ 



for 



ra 



pole 



140 GeV to ~ 1 for 



ra 



pole 



190 GeV. By increasing tan f3 one gradually departs 



from the |//| — > oo limit, and the enhancement to A^o from left-right mixing decreases. Thus, 



there is an interplay between A$ t and the overall factor of m\. (Of course, m h0 grows with 
tan (3 and hence with m p t ole .) 

As was pointed out in the previous section (see also Appendix B), // and the amount 
of mixing Y are subject to color and/or charge breaking (CCB) constraints. In this section 
points which are consistent with (13), i.e., which have no negative-energy CCB minima, are 
marked by a circle - O. Points which are inconsistent with (13) and therefore have either a 
local or global negative-energy CCB minimum are marked by a "diamond"- O. The filled 
diamonds correspond to a global color breaking minimum (as determined using numerical 
procedures), which is (most likely) unacceptable. The open diamonds represent CCB minima 
which are only local (i.e., above the standard model minimum) and are therefore safe. 
Constraint (13) is always evaluated at the t-scalar scale, independent of the scale chosen for 
minimization. The same is true for the numerical search of global CCB (GCCB) minima. 

Calculating m h o is reduced (in our case) to a great extent to calculating A^o, i.e., the 
deviation from what would be the case if the matter were supersymmetric at the weak scale 
[e.g., mi = mi 2 = m t in (14) and (15)]. Two standard ways to perform the calculation are 
the effective potential method (EPM) (e.g., see Ref. [45,44,46,36]) and the renormalization 
group method (RGM) (e.g., see Ref. [47,48,36]). The two methods correspond to a "run and 
diagonalize" and "diagonalize and run" algorithm, respectively. The algorithms and their 
implementations are discussed and compared in greater detail in Appendix C. Here we would 
like to demonstrate that the EP and RG methods do agree (within their domains of validity), 
i.e., the results are not sensitive to the way one treats the threshold corrections. In particular, 
both methods exhibit a large (10 — 20%) two-loop ambiguity from the presence of some two- 
loop residual terms in the calculations (next-to-leading and leading two-loop logarithms, as 
well as finite terms). We explicitly identify and discuss such terms in Appendix C. The 
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residual logarithmic two-loop terms can be described in terms of a scale ambiguity (i.e., 
the scale dependence of the one-loop calculation is of two-loop order). Indeed, comparing 
the two methods and studying that ambiguity allow us to estimate the two-loop correction. 
However, the uncertainty remains. 

In the EPM the ambiguity is explicitly related to a scale dependence, as one is required to 
specify a subtraction scale Q (i.e., halt running at Q and diagonalize mass matrices). Indeed, 
including AV one -' oop (Q) corrects for one-loop leading logarithms, and thus guarantees a 
scale independence of the calculation, but only to that order (and up to wave-function 
renormalization). It is suggestive to take Q = Mz } i.e., the scale at which the physical 
inputs are given. Another attractive choice is a scale at which AV two ~ loop (Q) is minimized. 
We will not attempt to find such a scale rigorously, but estimate Q 600 GeV. (One indeed 
expects such a Q to be in the vicinity of some average scalar-quark scale, as for A.V one ~ loop 
in Ref. [49].) We will use the comparison between the two methods to justify our choice. 
(Of course, there is no one choice that is suitable for all sets of parameters.) 

In the RGM, which sums the leading logarithms to all orders if one integrates numerically, 
the ambiguity is more apparent in terms of orders in perturbation theory (i.e., the order 
at which the summation is truncated in approximate iterative solutions). However, there is 
also an explicit scale dependence from the choice of decoupling scales for the heavy fields 
(i.e., decouple heavy eigenstates, calculate RGE's, and run to the next decoupling scale). 
We will decouple the heavy Higgs doublet H at m^o 2.1 TeV and other fields at lower 
thresholds 10 . When decoupling H we include appropriate matching conditions [48]. 

The EPM is more straightforward to implement (and is therefore more common in the 
literature) while the RGM requires a careful consideration of decoupling scales and matching 
conditions. 

The scatter plot Fig. 3 displays the results of an EPM calculation with Q = Mz- (When 
increasing Q the soft mass parameters and m t decrease so that the EPM m h o is diminished.) 
The scatter plot Fig. 4 displays the results of a RGM calculation (performed numerically). 
Comparing them, one observes a 10 — 20 GeV difference between the respective Higgs boson 
masses (for a given point in the parameter space). However, this is not due to the different 
choices of an algorithm, but to the different residual two-loop dependences in the two calcu- 
lations. To illustrate the maximum ambiguity we used the t-quark pole (running) mass in 
the calculations presented in Fig. 3 (Fig. 4). (The ~ 5% difference in the t-quark masses is 
formally of two-loop order.) In both plots one can observe a separation into a heavier and 
a lighter Higgs boson branches. This is due to the enhancement to ra^o in some cases from 
tL — Ir mixing, and is further discussed in the following section. 

To further compare the two methods, it is useful to consider the following exercise. 
The RGM calculation could be modified so that leading logarithms are not summed to all 
orders, e.g., one could use the one-loop leading-logarithm approximation formula of Ref. 
[47]. Alternatively, we could minimize the unsummed \n(m^/Q) leading (and/or next-to- 
leading) logarithms in the EPM by choosing Q raj. In either case the residual terms in 
the different calculations are qualitatively brought on the same footing. We carried out that 



10 We checked that the details of the choices for those intermediate scales introduce a ~ 1% 
uncertainty. 
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exercise numerically, and Fig. 5 (Fig. 6) compare the EPM with Q = Mz (Q = 600 GeV) 
calculation to the RGM leading logarithm approximation (numerical calculation). Once the 
residual two-loop dependence of the calculation is properly adjusted, the two algorithms 
agree well, in particular for Q 600 GeV. [All calculations presented in Fig. 5 (Fig. 6) use 
the t-quark pole (running) mass.] 

This all points at the most important source of uncertainty: residual two-loop (and 
wave-function) scale dependences and a ~ 5% a s dependence of m t . Two-loop terms fuel 
the ambiguity, which therefore cannot be removed in a satisfactory fashion within a one-loop 
calculation. (Even though we could identify some two-loop terms, the consistency of their 
removal is doubtful). Expanding the calculation to higher orders is not straightforward. We 
discuss two-loop calculations in Appendix C. 

Figs. 3 and 4 define the theoretical ambiguity in determining the loop induced m^o in 
this scenario. While the former sets the upper bound (for a given m pole ), the latter is 
a more accurate leading-logarithm calculation (and uses properly the running t-quark 
mass). However, within a one-loop calculation, and in particular when next-to-leading 
logarithms and some finite terms are comparable to higher-loop leading logarithms (which 
is the case here - see Appendix C), neither method provides a more consistent calculation 
then the other. As we argue in Appendix C, the former is also an upper bound on two- 
loop calculations, which most probably diminish m^o. We do not think that any stronger 
statements are justified at present. We continue to discuss the prediction and its upper 
bound in the next section. 

A different issue is that of the constraints on the t-scalar left-right mixing enhancement to 
A/jO. In the EPM [e.g., Eq. (14)] the enhancement is a straightforward result of diagonalizing 
the CP-even mass matrix. In the RGM one has to define the decoupling scale for the t- 
scalars. If raj ^> m^ 2 then one would decouple only t\ while t 2 loops could still contribute to 
A/jO. (In practice, the effect is accounted for in the RGM by a set of appropriate matching 
conditions.) However, the amount of mixing T (and thus the enhancement) is strongly 
constrained by the physical vacuum. 

Requiring that the SM extremum is a minimum (i.e., no negative squared masses cor- 
responding to physical degrees of freedom), and furthermore, is the global minimum (i.e., 
there is no deeper minimum that does not conserve color and/or charge) strongly constrains 
T. Constraints imposed by the latter may be evaded 11 if the non-SM global minimum is 
separated from the SM local minimum by a tunneling time greater than the age of the 
universe [50]. Constraint (13) (naively) attempts to eliminate all, i.e., local and global, 
negative-energy CCB minima in that region of parameter space and is too strong. The 
points marked in the scatter plots by only a an open "diamond" survive a more careful 
analysis, i.e., they correspond to only local minima 12 . We find, using numerical methods, 
that for m p t ole < 165 GeV (165 GeV < m pole ) roughly 80% (15%) of the points that fail to 



11 This claim, however, does not take into account finite-temperature effects. 

12 In fact, all points are found to have at least a negative-energy CCB local minimum of type 
M 2 = - see Appendix B. However, we find that in our region of the parameter space (13) is a 
useful tool to identify minima that compete with the standard model minimum. 
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satisfy constraint (13) correspond to (probably unacceptable) GCCB minima of the (one- 
loop improved) tree-level scalar potential. The enhancement is then partially washed away 
by requiring consistency with the physical vacuum. The possibilty of coexistence of the two 
vacua (with a tunneling time greater than the age of the universe) is, at the present stage of 
our calculation, a remaining ambiguity. We further discuss those issues in the next section 
and, in greater detail, in Appendix B. 

V. THE PREDICTION FOR THE HIGGS BOSON MASS AND ITS UPPER 

BOUND 

In this section we present and discuss our numerical results. In section IV we discussed 
various methods to calculate the mass (to one-loop). The effective potential method (EPM) 
with a subtraction scale Q = Mz (e.g., see Fig. 3) puts the weaker constraints on ra^o and is 
the conservative choice for the upper bound. In Fig. 7 we present monte-carlo distributions 
of m h o using that method for 155 < m p t ole < 165, 165 < m p t ole < 175, and 175 < m pole < 185 
GeV. The upper bound can be read off the figures, i.e., ra^o < 101, 106, 114 GeV for m pole < 
165, 175, 185 GeV, respectively. (Note that these bounds are deduced from our monte-carlo 
calculations, and are not rigorous.) Nevertheless, we would like to stress that typically ra^o 
lies below its upper bound. Indeed, one cannot rule out the possibility (if m pole is heavy 
enough) that m h o is just outside the range that LEPII may cover. However, this is not likely 
on both grounds, the m pole distribution [Eq. (6)] and the m h o (theoretical) distribution (as 
reflected in Fig. 7). 

The role of the large left-right mixing in the t-scalar sector can be seen in Figs. 7a and 
7b, which exhibit two peaks in the mass distribution, a "logarithmic" and a "mixing" peak. 
The former is due to the diagonal entries in the t-scalar mass matrix [Eq. (12)] while the 
latter, which is at larger mass, is due to the mixing enhancement [Ag t in Eq. (15)]. The 
interplay between the mixing and the m\ factor in the loop correction A^o [e.g., in Eq. (14)], 
and between A^o and the tree-level mass mf (that grows with tan f3 and, therefore, with 
m pole ) can be seen by comparing Figs. 7a-b with Fig. 7c. The mixing enhancement plays an 
important role for smaller values of m pole (A/3 — > 0). In fact, that interplay generates for 
a fixed set of soft parameters a local minimum in the ra^o distribution (not shown in the 
figures) near m pole 165 GeV. 

The left-right mixing enhancement is strongly constrained when requiring that the min- 
imum of the full scalar potential is the physical one, i.e., no fields other than H° = v& own 
and H® = v uv have non- vanishing expectation values (see section III and Appendix B). 
Constraint (13) eliminates potentially dangerous points in the parameter space (marked by 
a "diamond" in Figs. 3 and 4), some of which, however, correspond to only a local color- 
breaking minimum (marked by an open "diamond"). Those which correspond to global 
minima are marked by a filled "diamond". Fig. 8 is the same as Fig. 7 except that points 
which correspond to global color and/or charge breaking (GCCB) minima are omitted (the 
shaded areas have no negative-energy CCB minimum, global or local.) The mixing peak 
has diminished significantly but has not disappeared. Its presence is even stronger for larger 
values of m pole (see appendix B). (Some caution is required: the cosmological history of the 
universe may accommodate such a global minimum.) 
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A different issue is that of the large two-loop ambiguity in the calculation pointed out 
in section IV. For example, in Fig. 9 we show the same distributions as in Fig. 7, but in 
which the running (rather than the pole) t-quark mass was used to calculate ra^o. (The 
running mass is ~ 5% smaller.) Also, the distribution corresponding to the renormalization 
group method (RGM) calculation (e.g., see Fig. 4), which we showed 13 in Fig. 6 to roughly 
agree with the EPM calculation (with Q mi), is shown in Fig. 10 for 155 < m p t ole < 165 
and 165 < m p t ole < 175, respectively. Again, the t-quark running mass is used. Also shown 
are the distributions when GCCB points are omitted (shaded areas) to be compared with 
the total areas (shaded and unshaded) in Fig. 8. Note the concentration at lighter masses 
compared to previous figures, even though the upper limits are only slightly changed. It 
should be stressed that had we not included the proper RGM matching functions [48] (which 
would roughly correspond to setting the left-right mixings to zero by hand) then the RGM 
distribution would sharply peak at the lower end of the ra^o range. Thus, the importance 
of the matching functions is obvious. 

All, Figs. 7 - 10, are consistent within the theoretical ambiguity. It is suggestive that a 
consistent two-loop calculation will reduce the upper bounds given by the EPM calculation 
(we choose conservatively Fig. 7), but probably not to lower values than suggested by the 
RGM (Fig. 10). Further study is required before any stronger conclusions can be drawn. 
Even so, for m pole ^175 GeV ra^o [in case (1)] is within the reach of LEPII, and possibly 
of the Tevatron (see Ref. [51]). If m pole is lighter, as suggested by electroweak precision 
data, our histograms imply that m h o may still be relevant for Higgs searches at LEPI. (In 
particular, if two-loop corrections reduce ra^o as suggested by Fig. 10.) 

VI. CONCLUSIONS 

To summarize, we worked in a constrained framework assuming (up to high and low-scale 
threshold corrections) a\ = a 2 = a 3 = «g, hj, = h T} and universal initial conditions. Also, 
we explicitly included the correlation between the t-quark mass and the weak angle from 
electroweak precision data. Within that framework, only tan f3 ~ 1 and large tan f3 regions 
are allowed. In the large tan/? region (which is allowed only for a large m pole , unless mj 
is allowed to get large finite corrections from superpartner loops) m h o is bounded by the 
triviality bound, and is probably not interesting for Higgs searches at LEPII. 

In the tan (3 ~ 1 region the tree-level m^ nearly vanishes, but grows with m pole . The loop 
correction to the mass is enhanced by large — Ir mixing terms. The latter are diminished 
with an increasing m pole but, as m pole grows, the overall factor of m\ can compensate for 
that effect. The enhancement due to the large mixing is further constrained by requiring 
the physical vacuum to correspond to the global minimum of the full scalar potential. We 
find that for parameters corresponding to large mixing there is usually a color breaking 
minimum, but in many cases this is only a (presumably harmless) local minimum. We also 



The RGM calculation sums the leading logarithms to all orders in perturbation theory, while 
in the EPM calculation with Q ~ higher order terms are, in principle, minimized, and do not 
need to be summed. 
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pointed out that a two-loop ambiguity in the calculation of the Higgs boson mass is enhanced 
in our case and can be as large as 20 GeV. All these effects lead to the conservative bound 
m^o < 100 (110) GeV for m p t ole < 160 (175) GeV. (We briefly discuss two-loop calculations 
in Appendix C, where we conclude that such a calculation will most probably strengthen 
the upper bound.) The Higgs boson mass in the class of models characterized by tan f3 ~ 1 
is then almost certainly within the reach of LEPII (and possibly within the reach of the 
Tevatron [51]). Furthermore, the Higgs boson is typically lighter than the upper bound, and 
may be light enough to still be discovered at LEPI. 

Regarding the t-quark mass, one requires for SU{2) X U(l) breaking (in this scenario) 
m poie > Q e y 5 anc [ f ur ther requiring tan (3 > 1.1 (which is the case if the divergent limit 
is to be avoided, i.e., if the divergence is not stabilized by model-dependent finite terms) we 
obtain m p ole > 155 GeV. Eq. (6) and, in particular, the structure of the radiative corrections 
to the Z — >■ bb vertex, however, suggest that m pole is not much larger. (The latter can be 
evaded only if both the t-scalar and the Higgsino are light, which is not very likely in our 
case.) Also, we were able to understand the smallness of the parameter space allowed by 
hb = h T in terms of a custodial symmetry. The symmetry (theoretical) argument suggests 
m pole < 185 GeV (A^ < 1). In Appendix D we briefly examine extended Higgs sectors, and 
point out that the Higgs mass is a possible probe of such extensions if m pole (and/or tan (3) 
is known. 

In Appendix B we point out that the full (one-loop improved) scalar potential generically 
possesses color and/or charge breaking negative-energy local minima. Analytic constraints 
that are often found in the literature attempt to eliminate all such local minima (i.e., are 
too strong) but by construction apply to only certain types of minima while ignoring others 
(i.e., are too weak). We give a general classification of the local minima and are able to 
further identify [in case (1)] those that can potentially compete with the standard model 
minimum. 
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APPENDIX A: THE HIGGS SECTOR CUSTODIAL SYMMETRIES 

In this appendix we examine some symmetry aspects of the SGUT solutions, and in 
particular, of case (1). Some features of the latter, as well as the distinction between 
the two cases, can be better understood in terms of the custodial SU(2)l + r approximate 
symmetry of the SM lagrangian [26]; i.e., turning off hypercharge and flavor mixings, and 
if h t = hb = h } then one can rewrite the t and b Yukawa terms as an SU(2)l x SU(2)r 
invariant, i.e., 
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(where in the SM H 2 = it 2 H*) and SU(2) L X SU(2) R -> SU(2) L+R for v ± 0. However, 
/it 7^ hb and the different hypercharges of t° L and 6^ explicitly break the left-right symmetry, 
and therefore the residual custodial symmetry. 

In the MSSM, on the other hand, Hi is distinct from H 2 and if v& own 7^ v up SU(2)l x 
SU(2)r — > U(1)t 3L +t 3R [27]. In particular, the SU(2)l + r symmetry is preserved if f3 = ^ 
but is maximally broken if f3 = -|. Ignoring hypercharge, we can then define for illustration 
two quantities, 

hf — hf, 

Bm = Vdown , (A2a) 

m t - m b 

„ tan 8 — 1 , , , , 

B H = m b , A2b 

m t - m b 

that measure the SU(2)l + r symmetry breaking in the matter and Higgs sectors, respec- 
tively. If we map the tan (3 — m p t ole plane onto the Bm — Bh plane, we find that cases (1) and 
(2) approximately correspond to (1,0) and (0,1), respectively. Thus, while the constrained 
SO (10) solution realizes m t 7^ mj as an intrinsic property of the vacuum (i.e., v up ^> Vdown) 
and the grand-unified symmetry protects Bm ~ [29], the tan/? ~ 1 case realizes the 
SU(2)l + r (explicit) breaking as a scaling phenomenon (i.e., Bm oc h{ txed ) and the vacuum 
respects the symmetry (up to loop corrections). 

The symmetry that is spontaneously broken in case (2) is a O4 X O4 symmetry. Haber and 
Pomarol observed [28] that for m 3 — > [which is the situation in case (2)] there is no mixing 
between Hi and H 2 and the Higgs sector respects O4 X O4 (up to gauge-coupling corrections). 
The symmetry is broken to O3 X O3 for v& own 7^ v up 7^ and the six Goldstone bosons are the 
three SM Goldstone bosons, A , and H ± . The symmetry is explicitly broken for g 2 7^ (so 
m H + = Mw) and is not exact even when neglecting gauge couplings (i.e., m 3 7^ 0). Thus, A 
and H ± are massive pseudo-Goldstone bosons, i.e., m 2 H+ — M^y m 2 A o = C X m\. However, 

C ^ sin 2(3 

can be large, which is a manifestation of the fact that O4 X O4 — > O4 X O3 
for Vdown = 0. (The limit m 3 — > corresponds also to a U(l) Peccei-Quinn symmetry [13].) 

Our main interest in this work is, however, case (1), where from (9b) the Higgs-sector 
SU (2) custodial symmetry constrains 

+ P 2 = m H 2 + P 2 = _m 3- (A3) 

This is realized by taking the limit \fi\ — > 00, provided \B\ \fi\ which in our formalism 
comes out as a prediction. (Indeed, many SGUT models can easily realize a large B 0} 
for example, see Ref. [52].) Relation (A3) results in a (tree-level) massless eigenvalue in 
the neutral CP-even mass matrix. (Note that the SM Goldstone bosons are also massless 
Goldstone bosons of the broken SU{2) X SU{2).) 

The smallness of is protected by the custodial symmetry, i.e., 7^ is explained 
by loop corrections. The broken symmetry in the Yukawa sector induces at the loop level a 
small breaking in the Higgs sector (e.g., m 2 H 7^ mjj). The magnitude of the effect can be 
estimated by integrating J tan (3/ tan f3, i.e., 

A fj (M z ) = A fj (M G ) - JL£^(fc* _ hl)dln^- « A fj (M G ) + 0.6, (A4) 
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where we neglected h T , as well as threshold corrections, took h t 1, hj, 0, and assumed 
Ap < 1. Thus, even though the parameter space in case (1) is severely restricted, we can 
understand the restriction in terms of a symmetry. 



APPENDIX B: CONSTRAINTS ON COLOR AND CHARGE BREAKING 

MINIMA 

Let us consider the super-potential 

W = h t t L t R H° 2 + fiH^Hl (Bl) 

where we have performed an SU (2) rotation so that H 2 + has no vacuum expectation value. 
(We do not distinguish our notation for a superfield from that of its scalar component.) It 
is easy to convince oneself that only those terms in the super-potential are relevant when 
searching for the global minimum of the scalar potential in the large h t limit [53]. Using 
standard techniques we arrive at the corresponding scalar potential: 

h 2 V = M 2 — T + A, (B2) 

where the bilinear, trilinear, and quartic terms are 

M 2 = m\Ef + m\Bf + 2m 2 3 H°H° + m\ P L + m\ P R , (B3) 



r 



2 (\A t \H° + s\fi\H°) tj. 



(B4) 



anc 



A = (P L + t 2 R )H? + P L P R + 1(| x «D - terms- 



IBS) 



respectively. All fields were scaled <f> — > <f>/h t and are taken to be real and positive (our 
phase choice for the fields, which fixed m 2 < and T > 0, i.e., maximized the negative 
contributions to V) and all parameters are real. m\ 2 = m 2 H + /i 2 , s = fiA t /\fiA t \ } and the 
expression for the "D— terms" is [53] 

12 



— oli 
5 



'226 



2t 2 R 



+ a 2 [Hf-Hf + P L ] 2 + U [P L -P R ] 



(B6) 



We can parametrize the four fields in terms of an overall scale X and three angles < 
a, /3, 7 < -|: H° = X sin a cos /3, H 2 = X sin a sin /3, t R = X cos a cos 7, = X cos a sin 7, 
and redefine 



h 2 V(X) = M 2 (a, /3, 7 )X 2 - r(a, /3, 7 )X 3 + A(a, /3, 7 )X 4 . 



(B7) 



Then, for fixed angles, V(X) will have a minimum for X 7^ provided the condition 
32M 2 A < 9T 2 is satisfied. In that case, 
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x --- 

mm ~ 8 A 



1 / 32M 2 A N 



2 



> 0, (B8) 



and 



h]V mm = - l -Xl m {^-X mm — M 2 ) . (B9) 



We only allow parameters for which a standard model (SM) minimum exists. The SM 
minimum corresponds to a = | and f3 = (3° (7 is irrelevant and tan (3° = v up /vdown is the 
angle used to fix //, m 3 , as well as the Yukawa couplings). It is easy to convince oneself that 
in that limit the 4x4 second-derivative matrix is 2 X 2 block diagonal (otherwise baryon 
number is violated). Thus, it is sufficient to confirm that the four physical eigenvalues are 
positive to ensure that it is a minimum. (This is done in our numerical calculations.) If these 
conditions are satisfied then the SM is at least a local (negative-energy) minimum, and one 



has T SM = 0, Mj M < 0, and (B8) and (B9) reduce to the usual results X s ™ n = ^-M| M /2A SM , 

h*V*? n = -M| M /4A SM . 

Let us now consider the possibility of additional CCB minima with cos a 7^ 0. Since 
r > for the class of models we are considering this requires r / 0, which we assume 
hereafter. From (B7) - (B9) it is easy to classify the possible color and/or charge breaking 
(CCB) minima for definite a, f3, 7. One finds that for T 2 < 32AM 2 /9 there is no CCB 
minimum, while for 32AM 2 /9 < T 2 < 4AM 2 the CCB minimum exists but has V^°f > > 
V™ n , which is presumably safe. For 4AM 2 < T 2 (including the more rare case M 2 < 0, that 
must fall in this category) there is a negative- value CCB minimum, which may however be 
either local (presumably safe), i.e., V^™ > V^f n7 or global (probably unacceptable), i.e., 
^min < Kniri- Here, a sufficient (but not necessary) condition for an acceptable model is 

T 2 < 4AM 2 , (B10) 

which generalizes that of Ref. [53]. In principle, the above discussion holds for any number 
of fields (i.e., any number of angles), only the explicit expressions for M 2 , T, and A are 
more complicated. If constraint (B10) holds for every choice of a, f3, and 7 (cos a 7^ 0) 
then there is no negative- valued color and/or charge breaking minimum, global (GCCB) or 
local. In the special case M 2 < the constraint cannot be satisfied and there will be a 
negative-energy local CCB minimum, an observation made previously by Gunion et al. [53]. 

If (B10) does not hold, further investigation is needed to determine whether the minimum 
is global or local. It is straightforward to show that (B10) is satisfied for sin f3 = 0, i.e., there 
is no negative- valued CCB minimum for H® = for the class of models we are considering. 
We will therefore restrict our attention to the case H® 7^ 0, in which case it is convenient 
to reparametrize X = and rescale all fields by an additional factor of l/H®- This 
simplifies the discussion and enables us to examine and generalize previous work on the 
subject [54,43,55,53]. 

By inspection (B10) cannot hold for all values of the fields and there is always (and not 
just for M 2 < 0, which is hard to achieve in dangerous directions) a local CCB minimum, 
e.g., for M 2 = (A is positive definite). For example, one can choose 

H° = -4, (BHa) 

mi 
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t R = t L -^, (Bllc) 

m tR 

which is always possible because of (8). We confirmed this in our numerical studies, i.e., 
there is a CCB (negative- valued) local minimum even if V SM = V(vdowm v up) is the global 
minimum. (We do not consider additive constants.) Although such local minima may be of 
cosmological interest, the relevant question in our case is whether V SM is the global minimum 
of the full (or more precisely, of the four-field) scalar potential. Similarly, we assume that 
positive- valued CCB minima would not be populated and are therefore harmless. The answer 
has to be given by numerical mapping of all minima, but it is still useful to review how one 
could derive analytic constraints from (BIO), which are typically relevant only for specific 
regions of the parameter space. 

For example, if we fix H° = 0, = Ir = 1, then (BIO) gives the well known result [54] 

A 2 t <3(ml+ml + m 2 H2+l j 2 ). (B12) 

The equal field direction was chosen so the positive D— term contribution oc l/h 2 vanishes. 
This is not a relevant requirement for h t 1 (i.e., there may exist a deeper minimum with 
non vanishing D— terms). It is obvious that (B12) is not relevant when \fi\ — > oo. We find 
that a more useful constraint is [taking H° ~ 1, ~ Ir = t <C 1 and using (BIO)] 

(|A t | + 5 H) 2 <2(m? i+ my. (B13) 

Note that if the order t 2 corrections to (B13) are not negligible then V is less negative, which 
motivated our choice f < 1, i.e., the more dangerous direction. That constraint was also 
used by Drees et al. [43] to describe their numerical results (i.e., if we take m\ = m 2 in their 
formula). We confirmed in our numerical studies that (B13) and not (B12) is relevant in 
our case. [We found nearly all points to be consistent with (B12) for our ranges of the soft 
parameters.] 

Finally, we perform numerical minimization of the potential (B2) and map all minima 
using monte-carlo routines. We find that (B13) is to a good approximation sufficient to 
avoid a GCCB minimum, but is not necessary. For example, only about 80% (15%) of 
the points which are inconsistent with (B13) correspond to GCCB for m p t ole < 165 GeV 
(m p t ole > 165 GeV). We compare the points consistent with (B13) and with the numerical 
search in Fig. 8. The greater relevance of (B13) for m pole < 165 GeV can be understood if we 
recall that the lower bound on | cos 2/3| increases with m pole , and thus V SM becomes a deeper 
minimum. (The lower bound is related to h t unitarity, and that statement is independent 
of our hb = h T assumption.) On the other hand, \fi\ and the negative terms in V CCB [e.g., 
Eq. (B2)] diminish with increasing m pole . That interplay shifts the CCB minimum from a 
global to a local minimum. 

Let us demonstrate that claim in a simple-minded example assuming h t w h{ lxed w 1-1.1. 
We define sin 2 /3 = 2m 2 /h 2 v 2 = x, 0.5 < x < 1 \v = 2Mw 1 92\- We can rewrite 
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h 2 t V SM = - l -Mlhy [1 - 4x(l - x)] . (B14) 

|/z 2 V SM | sharply increases with m p t ole (i.e., with x). On the other hand, taking m 2 H ~ M\ sa 
in (9a) we have 

? » "'«, fe) • (BIS) 

and fj 2 decreases from +oo to in the above range ol x. Now assume that constraint (B13) 
is not satisfied, i.e., 

h*V™ « -±«V " ml - ml) t 2 , (B16) 

where t 4 <C ^ 2 and /c is a positive number in the range [|, |] (X m i n = k\[J,\). We also assumed 
\A t \ <C |//|, and neglected t 4 terms. (M 2 < would increase n.) To further simplify, consider 
the case m = 0, so that m~ ftim- sa 6M 2 sa 12m 2 rr . V CCB < V SM in this case if 

' l£ 'fl 2 

i ^i^ 2 ^ 2 < f i-M [« (a) - 24 1 fB17) 

4k 2 m^t 2 ~V2x-l7 l-4x(l-x)' 1 ; 

The r.h.s. has to be positive, i.e., x — > 0.5 + or m p t ole ^ 150 GeV. Slightly increasing m pole 
would reverse the inequality. 

Lastly, there are two caveats. The first is that we used in this section the (one-loop 
improved) tree-level potential. However, following Ref. [49] we perform the calculations at 
the t-scalar scale to eliminate large loop corrections. Secondly, a GCCB minimum may be 
"safe" if separated from the local standard-model minimum by a tunneling time greater than 
the age of the universe [50]. Such considerations are beyond the scope of our present work 
and would also require a consideration of finite-temperature effects. 



APPENDIX C: CALCULATION OF THE LOOP-INDUCED MASS 

1. The EPM: Run and diagonalize 

A straightforward way to calculate A^o is (a) to take the first and second derivatives of 
AV(Q) with respect to the neutral CP-even and CP-odd components of Hi at the minimum; 
(6) to absorb the latter in m^o; and (c) to calculate the correction to (11): the effective 
potential method (EPM). Like Ref. [44], we will follow Ref. [45], however, we add to their 
expressions "D-term" [46] and Higgs-Higgsino and gauge-gaugino [36] contributions, which 
are typically —(1 — 2) and —(2 — 4) GeV, respectively. Note that (z) all the parameters 
are taken at the subtraction scale Q. The mass matrices are thus calculated at the scale 
Q and (ii) only then are diagonalized, i.e., a "run and diagonalize" algorithm. Our two 
choices in section IV were Q = Mz and Q = 600 GeV ~ raj. (Large In logarithms 
are multiplied by small couplings and do not invalidate the loop expansion, i.e., the EPM. 
This is equivalent to the statement above that Higgs-Higgsino contributions are small.) One 
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corrects for threshold effects by the subtraction of one-loop leading logarithms (included in 
AV 7 ), and thus to that order the two choices are equivalent. 

The Q dependence of m^o^Q) comes about from (i) logarithmic wave-function de- 
pendence of Vi(Q) and (ii) two-loop implicit field-independent dependences of AV(Q) = 
AV[nii(Q), cti(Q), h t (Q)] (where ra 4 - runs over the relevant mass parameters). Both are of- 
ten underestimated. The importance of (i) is obvious, e.g., from the m\jM\ factor. The 
field-independent dependence (ii) was studied in Ref. [56], where it was shown to be of 
the order of two-loop (next-to-leading) logarithms. One can instead examine directly the 
expression for A^o, e.g., the leading logarithm in (14) evolves with scale (neglecting A$ t and 
multiplicative factors) roughly as 
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where k = 1,2 for Il,r, respectively, and we took m t (Q) ~ rn t (Q ), etc., raj ~ m t R 7 an d 
neglected all terms aside from h t and g 2 s (= 47ra s ) ones. (M g is the gluino mass.) This 
is of course a rough estimate, as, in practice all the parameters scale with Q. However, it 
illustrates qualitatively the presence of residual two-loop logarithms which are enhanced by 
the large couplings and masses when evaluating A^o. Note that typically the QCD term 
wins and the correction is positive near Q = Mz- 

Though formally of two-loop order, the residual Q dependence of the EPM ra^o is (in 
our case) of order 10 — 20%. The ambiguity in ra^o is reduced for mf ^ 0. Also, for values 
of tan f3 outside the allowed region A$ t is diminished (e.g., in the notation of Ref. [45] only 
A 2 2 is important), h t < 1 (which also diminishes the wave-function renormalization), and 
thus, there is a smaller ambiguity not only in ra^o but also in A^o. Nevertheless, one should 
be aware of such ambiguities, which are generically present in a calculation of that sort. 



2. The RGM: Diagonalize and run 

One could alternatively correct for thresholds by explicitly decoupling a particle below 
its threshold and redefining the effective field theory: the renormalization group method 
(RGM). Following Ref. [36,47], we diagonalize the Higgs mass matrices at scale Q = m A o } 
integrate out the heavy Higgs doublet, and define a one Higgs doublet model (1HDM) 
effective theory with a quartic coupling ^h 4 . Note that m^o ^> raj, so only the heavy Higgs 
doublet is integrated out at Q = m^o , and other heavy particles will be decoupled at lower 
thresholds. The boundary condition for X(m A o) is determined by the MSSM Higgs couplings 
at that scale and requires knowledge of tan /3(m A o) (i.e., wave-function renormalization). 
Mixing effects are accounted for by a complicated set of matching conditions [48,47]. A is 
then evolved down to Mz, where ra^o = \f\v = 2\/\Mw / 92', i-e., a "diagonalize and run" 
algorithm. 

One can solve the RGE's iteratively to get a semi-analytic approximation [36,47]. This is 
not an exact procedure, and at one-loop is accurate to leading logarithms (i.e., the accuracy 
of the EPM). Alternatively, if the RGM calculation is done by exact (or numerical) inte- 
gration, then the leading logarithms are summed to all orders. In agreement with Ref. [36] 
we find that X(Mz) (and thus ra^o) is diminished when integrating the RGE's numerically. 
However, the effect of the numerical integration (Fig. 5) is more important than anticipated 
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in Ref. [36] since A is heavier than the mass-scale used there and because of the large cou- 
plings. The difference between the numerical treatment and the approximate formula is of 
order of two-loop (leading) logarithms. The summation of the leading logarithms modifies 
the scale dependence of the calculated mass. 

A different scale dependence is via the details of the decoupling of heavy states. To ease 
the calculation, we do not diagonalize scalar-quark mass-squared matrices at an intermediate 
scale. Instead, we put in the appropriate matching conditions [48,47] by hand. These 
conditions are derived by expanding the EP assuming \m t u\ ^ m 2 ~ , etc., which still holds 
in our case [m f (m^o) < m t (Mz)\- Nevertheless, we repeated the IHDM RGM calculation 
but using the EPM with Q m A o to determine m h o(m A o) } and thus [using v(m A o)] the 
boundary condition for X(m A o). The agreement with the calculation using the boundary 
conditions of Ref. [48,47] is good. 

3. Validity of the comparison between the EPM and RGM 

Whether "running and diagonalizing" or vice versa leads to a slightly different answer. 
However, it was already noted in Ref. [36] that the difference between the two algorithms 
is small and of higher order, and can be ignored for the purpose of our discussion here. 
However, the inclusion of the appropriate RGM matching conditions is essential for the 
comparison between the two methods. 

More importantly, the RGM approximation of IHDM breaks down if is not sufficiently 
large. (In such cases the more conventional 2HDM RGM is appropriate.) Too small scalar- 
quark (diagonal) masses break down the expansion of Ref. [48,47] and too large scalar 
threshold-corrections can break down the EPM approximation. Indeed, one has to ensure 
that the scalar quarks are not heavier than about 1 TeV in order not to invalidate the loop 
expansion which underlies the EP, as these large logarithms would be multiplied by large 
couplings 14 (see, for example, Ref. [57]). We checked that the points in the parameter space 
where the RGM and EPM are in sharp disagreement (e.g., in Figs. 5 and 6) correspond to 
either one of these cases. 

4. Two-loop calculations 

In this work we do not carry out a two-loop calculation of A^o . Such a calculation is 
an elaborate task and attention should be paid to various issues, e.g., wave-function renor- 
malization and the coupled running of the masses and couplings; pole masses vs. running 
masses; the multi-scale structure of the low-energy theory; the correct choice of matching 
conditions (logarithmic and non-logarithmic) given a choice of decoupling scales; one-loop 
finite pieces (not included in the EP) that may become important in certain regions of the 
parameter space; in addition to correct counting of powers of the couplings. In fact, an 
hybrid algorithm between the EPM and RGM is required. 



This requirement is consistent with our choice of ranges for the soft parameters. 
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We received recently two interesting papers where (partial) two-loop analyses were pre- 
sented [58,59]. In this section we compare those results with the residual two-loop terms we 
identified above. We will take a somewhat critical point of view emphasizing that though the 
above studies (as well as the discussion here) pave the way to understanding A^o beyond 
one-loop order, they are not complete. We conclude that the one-loop leading-logarithm 
result is most probably an upper bound on the two-loop result, but we doubt whether any 
stronger conclusions can be drawn at the present. 

Above, we presented the leading-logarithm result. In the EPM this is given by /\y one - lo °p 
and in the RGM one solves the RGE's to that order. We found no significant difference 
between the two methods. Next, we proceeded to calculate beyond the leading order and 
integrated the RGM RGE's numerically, summing the leading logarithms to all orders in 
perturbation theory, and, in particular, to two-loop order. We then found a ~ —10 — 20% 
correction to A^o. On the other hand, we identified in the EPM expressions next-to-leading 
two-loop logarithmic terms. Those terms come about from the scale dependence of the 
running parameters, and introduce a similar (but positive) correction to A^o. Also, we 
noted the important (~ 10%) effect of m\ ole / m r t unnm3 [I + (4a s /37r)] [e.g., ( m ™ nnmfl ) 4 ~ 
0.8(mf°' e ) 4 ]. It is then suggestive that the two-loop leading logarithms and m t corrections 
would diminish the one-loop (leading-logarithm) result, while the next-to-leading logarithms 
would partially counter-balance that effect. Thus, we reach the conclusion that the one-loop 
leading-logarithm prediction for ra^o (which contain some next-to-leading logarithms) is an 
upper bound on the two-loop prediction. Any stronger conclusion would require a detailed 
consideration of the issues listed above. Furthermore, the significance of a detailed two- 
loop study in the context of SGUT's is not clear unless complemented by similar studies of 
threshold effects near the high-scale boundary. 

Let us proceed with the issue of the orders in perturbation theory. It was already pointed 
out by Chankowski [57] that the inclusion of /\y one - lo °p [ n (7) with only one-loop RGE's for 
the (RG-improved) potential parameters is an inconsistent procedure when counting powers 
of couplings. These issues have recently received great attention in the context of theories 
with a single scalar field (i.e., </> 4 or the SM) [60,61]. It was explicitly shown that a RG- 
improved L-loop effective potential is consistent when the parameters are calculated using 
(L+l)-loop RGE's [61]. 

The observations of Ref. [60,61] were recently applied to the MSSM by Kodaira et al. 
[58] who extended a previous work of Espinosa and Quiros [62]. Indeed, they find that the 
two-loop RGM with one-loop EP slightly increase the one-loop RGM result. However, they 
impose the naive assumption of only one relevant scale (Msusy) below which the SM RGE's 
are in effect (and the EP contain only a m t correction). Thus, the MSSM parameters are 
not treated on the same footing as the SM ones, and one expects some modification of their 
numerical results. However, we note that the correction is positive: recall that the difference 
between the two RGM approximations is the next-to-leading logarithms, which we found to 
be positive by examining the scale dependence of the result. 

Hempfling and Hoang [59] explicitly calculate the two-loop RGM correction (using the 
methods of Ref. [47]), finding 
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which roughly agrees with our Eq. (CI). 
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Those authors go further to derive a diagramatic two-loop result (in a simplified model) 
which diminishes A^o below its two-loop RGM prediction. They then show that, e.g., if 
the difference between the pole and running t-quark mass is accounted for in the RGM 
calculation by appropriate matching conditions, the two methods agree. This again stresses 
the need for adequate choice of matching conditions. However, it is not clear that there are 
no other effects of that order of magnitude that can slightly increase A^o. [Also, we showed 
that scalar mixings which are neglected in the above two-loop calculations and in our Eq. 
(Cf), are not apriori negligible.] 

In short, all analyses support the conclusion that two-loop logarithms will diminish the 
one-loop leading logarithm result (in particular, if one properly distinguishes the t-quark 
pole mass from the running mass). However, it is not clear that the prediction of a complete 
two-loop analysis lies as low as suggested in Ref. [59] (or by our Fig. 10). (If indeed this is 
the case our upper bounds will be significantly strengthened.) 

APPENDIX D: EXTENDED HIGGS SECTORS 

If there is an extended heavy Higgs sector (as is assumed by some authors when trying to 
explain the light fermion spectrum) with large representations coupling to the third family 
(i.e., there is no flavor symmetry that forbids such couplings) then (2) does not hold and our 
constraints are evaded. Alternatively, some models extend the light Higgs sector by adding a 
SM singlet superfield, S, to the spectrum. This is the only addition to the Higgs sector which 
is consistent (i.e., without fine-tuned cancellations) with coupling constant unification. 

The new singlet (s) would mix with the other neutral Higgs bosons and would thus con- 
tribute to mfo', e.g., [63], 

9 v 2 
mlo < M| cos 2 2/3 + A 2 y sin 2 2/3, (DI) 

where A s (/c s ) is a new Yukawa coupling in the superpotential: X s HiH 2 S (k s S 3 ). (We do 
not distinguish here our notation of superfields from that of their scalar components.) The 
MSSM limit with ?tj^o — is recovered if A s — > 0, which is indeed the case if h t « h{ lxed [64]. 
However, the new Yukawa couplings (i.e., A s 7^ 0) enable one to realize (2) for somewhat 
smaller values of h t < h{ txed (and tan/? > I). The shift in the allowed region is shown in 
Fig. II for A s = 0.5 and k s = 0. Thus, mf can be large and ra^o can be heavy ~ 150 GeV 
[19,20,64] when loop corrections are added. (See also Ref. [65].) 
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FIG. 1. The i-quark pole mass (in GeV) - tan/3 plane is divided into five different regions. 
Two areas (low- and high-tan /3 branches) are consistent with perturbative two- Yukawa (b and r) 
unification and with an upper bound of ^|| GeV for the 6-quark current mass. The 0.85 factor takes 
into account a ~ ±15% theoretical uncertainty in the numerical calculation due to the matching 
conditions. Between the two branches the 6-quark mass is too high. For a too low (high) tan/3, 
the i-quark (6-quark) Yukawa coupling diverges. The region where all three (third-family) Yukawa 
couplings unify (dashed line) intersects the allowed high-tan /3 branch. [The three- Yukawa region 
is calculated without imposing any constraints on mj, and is assigned a ~ 5% uncertainty (shaded 
area) from corrections to the hf/hf, ratio.] The i-quark mass range suggested by the electroweak 
data and the tan /3 = 1 line are indicated (dotted lines) for comparison. 
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FIG. 2. The prediction for the Higgsino mass parameter ^ (after absorbing AV in the redefined 
tree-level parameters) as a function of (a) the i-quark pole mass and of (b) tan/3. (The i-quark 
mass and tan/3 are strongly correlated.) Filled squares indicate (probably unacceptable) points for 
which the SM minimum is only a local minimum. All masses are in GeV 
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FIG. 3. A scatter plot of the Higgs boson mass (calculated in the EPM with the subtraction 
scale Q at the Z-pole) vs. the i-quark pole mass. The o indicate points which have no (potentially 
dangerous) negative-energy color and/or charge breaking minimum. The o have such a minimum 
but it is a (safe) local one, i.e., it lies above the standard-model minimum. Points indicated by a 
filled diamond have an unacceptable CCB global minimum. The Z mass is indicated for comparison 
(dashed line). All masses are in GeV. 
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FIG. 5. The EPM (with the subtraction scale Q at the Z-pole) calculation of the Higgs boson 
mass is compared with the mass calculated in the RGM but using the Haber and Hempfling leading 
logarithm formula. All masses are in GeV. Points corresponding to global GCCB minima (indicated 
by a filled diamond in Figs. 3 and 4) are omitted. 
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FIG. 6. Same as in Fig. 5, except Q = 600 GeV and the RGM calculation is carried out 
numerically. 
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(c) 175 <m[ < 185 GeV 
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FIG. 7. The Higgs boson mass (in GeV) distribution in a sample monte-carlo calculation (using 
the EPM with the subtraction scale Q at the Z-pole) for the i-quark pole mass in the range (a) 
[155, 165], (b) [165, 175], (c) [175, 185] GeV. 
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FIG. 8. Same as in Fig. 7 except for the omission of points which correspond to a global color 
and/or charge breaking minimum (marked by a filled diamond in Fig. 3). For comparison, the 
shaded areas indicate points which are consistent with the constraint (13) (marked by a circle in 
Fig. 3). The difference between the two distributions is points with only a local CCB minimum 
unnecessarily excluded by (13) (marked by open diamonds in Fig. 3). 
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FIG. 9. Same as Fig. 7 except for substituting the i-quark running (rather than pole) mass in 
the mass expressions for the scalars. 
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FIG. 10. Same as in Fig. 7a-b except using the RGM (with numerical integration). The dis- 
tributions omitting points which correspond to a global color and/or charge breaking minimum 
(marked by a filled diamond in Fig. 4) (lightly shaded areas) are also indicated, and should be com- 
pared with the total (shaded and unshaded) areas in Fig. 8. (For a i-quark heavier than about 180 
GeV the 1HDM assumption that underlies the calculation is not always accurate - see Appendix 
C.) 
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FIG. 11. The low tan/3 branch in the presence of an additional Higgs field with A s = 0.5 and 
k s = 0. The allowed area when there is no singlet field (A s = 0) is shown for comparison (shaded 
region), the i-quark pole mass (given in GeV) range suggested by precision data and tan/3 = 1 are 
indicated (dotted lines). 
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